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Abstract. We have investigated the infrared spectra of all 46 Herbig Ae/Be stars for which spectroscopic data are available 
in the ISO data archive. Our quantitative analysis of these spectra focuses on the emission bands at 3.3, 6.2, "7.7", 8.6 and 
11.2 micron, linked to polycyclic aromatic hydrocarbons (PAHs), the nanodiamond-related features at 3.4 and 3.5 micron, the 
amorphous 10 micron silicate band and the crystalline silicate band at 11.3 micron. We have detected PAH emission in 57% 
of the Herbig stars in our sample. Although for most of these sources the PAH spectra are similar, there are clear examples 
of differences in the PAH spectra within our sample which can be explained by differences in PAH size, chemistry and/or 
ionization. Amorphous silicate emission was detected in the spectra of 52% of the sample stars, amorphous silicate absorption 
in 13%. We have detected crystalline silicate emission in 1 1 stars (24% of our sample), of which four (9%) also display strong 
PAH emission. We have classified the sample sources according to the strength of their mid-IR energy distribution. The systems 
with stronger mid-infared (20-100 yum) excesses relative to their near-infrared (1-5 yum) excess display significantly more PAH 
emission than those with weaker mid-infrared excesses. There are no pronounced differences in the behavi our of the silic ate 
feature between the two groups. This provides strong observational support for the disk models by Dullemo nd et alJ yOOl|), in 
which systems with a flaring disk geometry display a strong mid-infrared excess, whereas those with disks that are strongly 
shadowed by the puffed-up inner rim of the disk only display modest amounts of mid-infrared emission. Since the silicates are 
expected to be produced mainly in the warm inner disk regions, no large differences in silicate behaviour are expected between 
the two groups. In contrast to this, the PAH emission is expected to be produced mainly in the part of the disk atmosphere that 
is directly exposed to radiation from the central star. In this model, self-shadowed disks should display weaker PAH emission 
than flared disks, consistent with our observations. 

Key words, circumstellar matter — stars: pre-main-sequence — planetary systems: protoplanetary disks 



1 . Introduction ISO played a major role in opening up the field of infrared 

j^ spectroscopy. The emission features at 3.3, 6.2, 7.7, 8.6 and 

■ - ' Although the optical— sub-mm energy distribution of Herbig 112 micron, found in ISO spectra of many HAEBE stars and 

Ae/Be (HAEBE) stars has been well explored by previous au- previously known as the unidentified infrared (UIR) bands, 

thors (e.g. Hillenbrand et al. 1992), the chemical and miner- ^^ generally attributed topolycyclic aromatic hydrocarbons 

alogical composition of the dust remained poorly studied un- (p^Hs, Leaer & Pug'^ ll984t Ullamandola et al. 1989). PAHs 

til the 1995 launch of the Infrared Space Observatory (ISO, are a large family of molecules, for which the fundamental 

[Kessler et al. 1996). This first possibility to study the com- ingredients are polycycHc benzene rings. PAH molecules are 

plete infrared spectrum of these objects in detail revealed a thought to be excited by far-ultraviolet (UV) photons. The ab- 

large variety in dust properties, from small aromatic hydro- sorption of such a photon induces a transition of the PAH 

carbons to siHcate dust. Moreover, some sources were shown molecule to an upper electronic state. The excited molecule 

to contain partially crystalline dust grains, similar to those then makes rapid transitions to a lower electronic state, leav- 

found in comets in our own solar system ("Waelkens et al.' ing most of the initially absorbed energy in the form of vibra- 

.1996; Malfait et al. 1998. 1999; van den Ancker et al. 2000a b; tional energy in the CC and CH bonds. PAHs cool down by in- 

Meeus et al. | 2001| ). frared (IR) emission in bands linked to these vibrational modes 



(Peeters 2002, and references therein). The 3.3 pm band is due 

* Based on observations with ISO, an ESA project with instruments to CH bond stretching vibrations, the 6.2 pm band is Hnked 

funded by ESA Member States (especially the PI countries: France, to the CC stretching mode and the 7.7 pm band corresponds 

Germany, the Netherlands and the United Kingdom) and with the par- to a combination of the CC stretching and the CH in-plane 

ticipation of ISAS and NASA. 
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bending modes. The 8.6 and 11.2 jjm features are linked to 
the CH in-plane and out-of -plane bending modes respectively 
dVermeii et alJ2002l and references therein). 

In a few sources, emission bands at 3.4 and 3.5 micron 
are observed. These spectral features are attributed to the 
CH stretching modes o f hydrogenated nanodiam onds (NANs, 
iGuillois et aU1999.:.Van Kerckhoven et al.l2002h . 

The infrared spectra of HAEBEs also contain silicate fea- 
tures. The broad feature at 10 micron has been attributed to 
the SiO bonds in warm, small amorphous silicate grains like 
oUvine. The feature can appear in emission or in absorp- 
tion in HAEBE spectra. The spectral signature at 11.3 mi- 
cron, as well as many features at longer wavelengths which 
are not included in the analysis presented in this paper, 
is hnked to cryst alline silicates (Malfait et al. 1998, 1999; 
Ivan den Ancker et al-2000a.b:.Bouwman et al-200 0. 2001). 

The spectral energy distribution (SED) of HAEBE stars 
is characterized by the presence of an IR flux excess, due to 
thermal emission of circumstellar matter. The geometry of 
the circumstel lar dust has been t he sub j ect of a long - lasting 
debate (e.g. iHill enbrand et al.' 1992; BerrilH et al. 1992; 
Hartmann et alJ Il993: .Bohm&Catala .1994: .Grininetal 



1994 Il996t bradvetalJ 119951: ICorcoran & Ravi Il998t 
Mannings & Sargennll997l) . For late-B, A and F stars, the 
evidence for the presence of disk-like geometries is generally 
accepted. For early-B stars, the matter is less clear. In these 
systems, the dissipation time scale of the circumstellar spheri- 
cal envelope is of the order of the pre-main-sequence life time. 
Disks as well as spherical envelopes might be present. 

In their study of the ISO spectra of 14 isolated Herbig 
Ae/Be stars, Meeus et alJ ll200L henceforth MOl) classified 
their sample into two groups, based on the shape of the SED. 
Group I contains the sources in which a rising mid-IR (20- 
100/im) flux excess is observed; these sources have an SED 
that can be fitted with a power-law and a black-body contin- 
uum. Group II sources have a more modest mid-IR excess; 
their SEDs can be reconstructed by a power-law only. MOl 
suggest phenomenologically that this classification represents 
different geometries of the circumstellar disk: group I sources 
ha ve flared disk s, gr oup II members have flat disks. 

iDullemondl ll2002l henceforth D02) and 

iDuUemond & DominikI ll2004 henceforth DD04) have 
modelled young stellar disks with a self-consistent model 
based on 2-D radiative transfer coupled to the equation of ver- 
tical hydrostatics. The model consists of a disk with an inner 
hole (~0.5 AU), a puffed-up inner rim a nd an outer part. The 
outer part of the disk can be flared (as in lChiang & GoldreichI 
I1997 I). but can also lie entirely in the shade of the inner rim. 
The SEDs of flared disks display a strong mid-IR flux excess, 
while self-shadowed disks have a much more modest mid-IR 
excess. D02 explains quantitatively the difference in SED 
shape in HAEBEs (as expressed by the classification of MOl) 
as the result of a different disk geometry; group I sources have 
flared disks, group II sources have flat self-shadowed disks. 

In this paper we investigate a possible link between the 
shape of the SED — a proxy for the geometrical distribution 
of the circumstellar matter — and the strength and profiles of 
circumstellar infrared emission bands in HAEBE systems. All 



spectra of HAEBE stars taken by ISO are investigated, amongst 
others spectra that have never been published before. It is there- 
fore the most complete sample of near-IR spectra of HAEBE 
stars ever investigated as a whole. Our study is not only qual- 
itative, but contains also, as opposed to many other articles, a 
quantitative analysis of the infrared spectra of HAEBEs. 

2. The data set 

2.1. The sample 

Our h st of HAEBEs is based on the catalogue of The et alJ 
0994 henceforth T94). We restricted ourselves to Table 1 and 
2 of that article. We enl arged the sample with HAEBEs studied 
by iMalfait et aP (Il998l henceforth M98), that satisfy the cri- 
teria postulated by T94. The resulting list was cross-correlated 
with the ISO data archive; we selected all objects for which ISO 
spectra were available with an ISO pointing within a distance 
of ±5 arcseconds of the position given by T94 and M98. 

We eliminated some of the remaining sources from the 
sample; TY CrA was omitted be cause the infrared emissio n 
emanates from the "TY CrA bar" dSiebenmorgen et alJbOOOl) . 
which is not directly associated to the star. The object also 
does not display Ha emission, which indicates that this star is 
no longer actively accreting. From a comparison between the 
ISO-SWS spectrum of MWC 137 and pho tometric measure - 
ments from the MSX Point Source Catalog lE ganetalJll997t) . 
it appears that there is an offset of 12" between the position 
indicated in T94 and the position of the infrared source. Only 
a small portion of the object's flux is seen in the SWS spec- 
trum. We therefore discarded this spectrum from our analy- 
sis. A similar problem was encountered for the SWS spectrum 
of IRAS 12496-7650. A mispointing of the ISO instrument of 
about 17" compar ed to the 2MASS point-source-catalogue po- 
sition dCutri et al.l l2003 ) caused an anomaly between the ISO 
S WS-spectra and the IR photometry. However, the PHT-S spec- 
tra for this source were not affected due to the larger aper- 
fiire of the ISO-PHOT instrument (see Sect. 122). 51 Oph and 
MWC 300 were also removed fror n the list; these are probably 
evolved objects and not HAEBEs dvan den Ancker et al.l200ll: 
JMolster et al. 2002, respectively). 

The resulting sample of 46 sources is tabulated in TableQ] 
The columns contain respectively the object's name, the type 
of ISO spectrum (see Sect. 12. 2t . the ISO Observation Sequence 
Number (OSN), the right ascension (RA), the declination (Dec) 
(both epoch 2000 coordinates), the date and starting time (UT 
start) of the measurement and the total integration time T in 
seconds. The objects are sorted by increasing RA. The given 
coordinates are the pointing coordinates of ISO for that source. 

2.2. Data Reduction 

We retrieved all spectra of HAEBE stars present in the ISO data 
archiv e'. Three types of ISO d ata were use d in this study: ISO- 
SWS Me Graauw et al."l996'i, ISO-PHT JLemke et alJll996l) 
and ISO-CAM (Cesarsky et al. 1993)- 



http://www.iso.vilspa.esa.es/ida/index.html 
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Table 1. The sample of HAEBEs used in this study, based on lThe et al.l(ll994 and lMalfait et alJ ( fl99 
observed ISO spectra ai-e given; SOI refers to SWS AOTOl, P40 to PHT-S and C04 to CAM04 spectra. 



For each object, the 



Sample stars 



Object 



AOT 



OSN 



RA (2000) 
h m s 



Dec (2000) 



Date 



UT start 
h m s 



T 

[s] 



V376 Cas 


SOI 


43501514 


00 1126.6 


-H58 50 04 


24-Ian-1997 


17:00:14 


3554 


VXCas 


P40 


58704023 


00 3130.5 


-1-6158 51 


25-Iun-1997 


15:42:53 


364 


Elias 3-1 


SOI 


67301306 


04 18 40.7 


-H28 19 16 


19-Sep-1997 


01:03:47 


3454 


AB Aur 


SOI 


68001206 


04 55 45.7 


-1-30 33 06 


26-Sep-1997 


05:26:35 


3454 


HD 31648 


SOI 


83501201 


04 58 46.1 


-H29 50 38 


27-Feb-1998 


19:47:48 


3454 


UXOri 


P40 


85801453 


05 04 30.0 


-03 47 14 


22-Mar-1998 


19:37:29 


364 


HD 34282 


SOI 


83301240 


05 16 00.5 


-09 48 34 


25-Feb-1998 


18:57:46 


1912 


HD 34700 


SOI 


66302638 


05 19 41.4 


-1-05 38 43 


09-Sep-1997 


11:10:28 


1912 




P40 


63602294 


05 19 41.4 


-1-05 38 42 


13-Aug-1997 


12:51:00 


172 


HD 35187 


SOI 


69501139 


05 24 01.2 


+24 57 36 


lO-Oct-1997 


23:44:54 


1912 


BFOri 


P40 


70101958 


05 37 13.3 


-06 35 01 


17-Oct-1997 


08:33:28 


364 


RRTau 


P40 


86603163 


05 39 30.5 


-1-26 22 26 


30-Mar-1998 


21:38:08 


364 


ZCMa 


SOI 


72201607 


07 03 43.2 


-1133 07 


07-NOV-1997 


04:42:28 


3454 


HD 95881 


SOI 


10400818 


110157.8 


-71 30 52 


29-Feb-1996 


05:04:40 


3462 




P40 


10400919 


110157.8 


-7130 52 


29-Feb-1996 


06:03:04 


140 


HD 97048 


SOI 


14101343 


1108 04.6 


-77 39 17 


06-Apr-1996 


10:15:08 


3462 




SOI 


61801318 


1 1 08 04.6 


-77 39 17 


26-M-1997 


08:30:31 


6538 




P40 


07900309 


1 1 08 04.6 


-77 39 17 


04-Feb-1996 


06:06:00 


172 




P40 


14101580 


1108 04.6 


-77 39 17 


06-Apr-1996 


11:26:28 


364 




P40 


62501510 


1108 04.6 


-77 39 17 


02-Aug-1997 


04:36:23 


620 




C04 


71801836 


1108 03.5 


-77 39 17 


03-NOV-1997 


02:54:09 


2208 




C04 


71901688 


1 1 08 04.6 


-77 39 17 


04-NOV-1997 


04:17:02 


2846 


HD 100453 


P40 


26000131 


11 33 05.6 


-54 19 29 


02-Aug-1996 


17:38:16 


140 




SOI 


26000230 


11 33 05.7 


-54 19 29 


02-Aug-1996 


17:41:20 


1912 


HD 100546 


SOI 


07200660 


11 33 25.5 


-70 1142 


28-Ian-1996 


13:35:29 


1044 




SOI 


27601036 


11 33 25.3 


-70 1142 


07-M-1996 


11:37:28 


1912 




P40 


10400537 


11 33 25.7 


-70 1142 


29-Feb-1996 


04:14:02 


140 


HD 104237 


SOI 


10400424 


12 00 06.0 


-78 11 34 


29-Feb-1996 


03:40:48 


1834 




SOI 


23300524 


12 00 05.1 


-78 11 34 


07-Iul-1996 


11:37:28 


1912 




P40 


23300625 


12 00 05.1 


-78 11 34 


07-M-1996 


12:10:04 


140 




P40 


53300118 


12 00 05.1 


-78 11 34 


02-May-1997 


02:32:49 


364 


IRAS 12496-7650 


SOI 


23300112 


12 53 15.9 


-77 07 02 


06-M-1996 


19:27:51 


3455 




P40 


07901717 


12 53 16.1 


-77 07 02 


04-Feb-1996 


12:56:06 


140 


HD 135344 


SOI 


10401575 


15 15 48.4 


-37 09 16 


29-Feb-1996 


10:10:12 


1834 




P40 


10401742 


15 15 48.4 


-37 09 16 


29-Feb-1996 


11:47:12 


140 




P40 


10401876 


15 15 48.4 


-37 09 16 


29-Feb-1996 


12:11:18 


140 


HD 139614 


SOI 


29701542 


15 40 46.3 


-42 29 53 


09-Sep-1996 


03:28:45 


1912 




P40 


10402322 


15 40 46.5 


-42 29 55 


29-Feb-1996 


14:45:38 


140 


HD 141569 


SOI 


62802937 


15 49 57.6 


-03 55 16 


05-Aug-1997 


04:40:44 


1912 




P40 


62701662 


15 49 57.7 


-03 55 17 


04-Aug-1997 


08:02:39 


236 


HD 142527 


SOI 


10402046 


15 56 42.1 


-42 19 24 


29-Feb-1996 


12:42:16 


1834 




P40 


10402547 


15 56 42.1 


-42 19 24 


29-Feb-1996 


15:42:20 


140 


HD 142666 


SOI 


10402952 


15 56 40.1 


-22 0141 


29-Feb-1996 


18:23:12 


3462 




SOI 


44901283 


15 56 40.1 


-22 01 39 


07-Feb-1997 


19:12:46 


6538 




P40 


10402847 


15 56 40.1 


-22 0141 


29-Feb-1996 


17:59:06 


140 


HD 144432 


SOI 


45000284 


16 06 58.0 


-27 43 08 


08-Feb-1997 


08:28:30 


6539 




P40 


10402662 


16 06 58.0 


-27 43 10 


29-Feb-1996 


16:01:26 


140 


HR 5999 


SOI 


28901506 


16 08 34.3 


-39 06 19 


01 -Sep- 1996 


02:46:48 


1912 




P40 


28901748 


16 08 34.2 


-39 06 18 


01 -Sep- 1996 


04:28:10 


140 




C04 


45800955 


16 08 34.3 


-39 06 19 


16-Feb-1997 


17:01:44 


2198 


Wra 15-1484 


SOI 


29901001 


16 27 14.9 


-48 39 27 


lO-Sep-1996 


20:04:23 


1912 


HD 150193 


SOI 


08200444 


16 40 17.9 


-23 53 45 


07-Feb-1996 


07:02:19 


1044 




P40 


64102335 


16 40 17.9 


-23 53 44 


18-Aug-1997 


13:13:35 


236 


AKSco 


SOI 


28902101 


16 54 44.8 


-36 53 18 


Ol-Sep-1996 


07:24:34 


1140 




P40 


64402829 


16 54 44.8 


-36 53 17 


21-Aug-1997 


13:00:15 


236 
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Sample stars (continued) 


Object 


AOT 


OSN 


RA (2000) 


Dec (2000) 


Date 


UT start 


T 








h m s 


o / " 




h m s 


[s] 


CD-42° 11721 


SOI 


08402527 


16 59 06.8 


-42 42 08 


09-Feb-1996 


16:26:56 


1816 




SOI 


64701904 


16 59 05.8 


-42 42 15 


24-Aug-1997 


11:11:16 


1912 




P40 


28900460 


16 59 06.8 


-42 42 08 


31- Aug- 1996 


17:58:54 


364 


HD 163296 


SOI 


32901191 


17 56 21.4 


-21 57 20 


lO-Oct-1996 


21:05:00 


3454 




P40 


32901192 


17 56 21.3 


-21 57 20 


lO-Oct-1996 


22:03:18 


140 


HD 169142 


SOI 


13601359 


18 24 30.0 


-29 46 50 


Ol-Apr-1996 


09:46:06 


1834 




P40 


13601437 


18 24 29.8 


-29 46 50 


Ol-Apr-1996 


10:30:40 


140 


MWC 297 


SOI 


70800234 


18 27 39.5 


-03 49 52 


23-Oct-1997 


18:56:23 


1912 


VVSer 


P40 


47800913 


18 28 47.9 


+00 08 40 


08-Mar-1997 


16:27:17 


364 


RCrA 


SOI 


14100458 


19 01 53.9 


-36 57 10 


06-Apr-1996 


04:48:02 


1834 




SOI 


70400558 


19 0153.9 


-36 57 10 


19-Oct-1997 


20:12:45 


1912 




P40 


11501230 


19 01 53.4 


-36 57 04 


11 -Mar- 1996 


11:04:57 


140 


TCrA 


SOI 


33402096 


19 01 58.8 


-36 57 49 


16-Oct-1996 


04:15:14 


3454 




SOI 


68900196 


19 01 58.8 


-36 57 49 


04-Oct-1997 


18:11:15 


3454 




P40 


14100562 


19 0158.8 


-36 57 49 


06-Apr-1996 


06:07:04 


140 


HD 179218 


SOI 


32301321 


19 11 11.2 


+ 15 47 17 


05-Oct-1996 


03:27:54 


3454 




C04 


36401031 


19 11 11.3 


+ 15 47 16 


14-NOV-1996 


16:11:58 


2240 


WWVul 


SOI 


17600305 


19 25 59.0 


+21 12 30 


lO-May-1996 


23:54:57 


1834 




P40 


17600465 


19 25 59.0 


+21 12 31 


11 -May- 1996 


01:33:17 


140 




P40 


51300108 


19 25 58.6 


+21 12 31 


12-Apr-1997 


03:50:24 


364 


BD+40°4124 


SOI 


35500693 


20 20 28.3 


+412151 


05-NOV-1996 


18:19:26 


3454 




P40 


15900568 


20 20 28.3 


+412151 


24-Apr-1996 


03:39:08 


140 


LkHa 224 


SOI 


85800502 


20 20 29.2 


+41 21 27 


22-Mar-1998 


09:49:31 


1912 


LkHa 225 


SOI 


85800403 


20 20 30.4 


+41 21 27 


22-Mar-1998 


08:51:15 


3454 


PVCep 


SOI 


14302273 


20 45 54.0 


+67 57 36 


08-Apr-1996 


13:37:16 


1834 




P40 


14302274 


20 45 54.0 


+67 57 36 


08-Apr-1996 


14:08:34 


140 


HD 200775 


SOI 


33901897 


21 01 36.8 


+68 09 49 


21 -Oct- 1996 


03:36:10 


3454 




C04 


10702605 


2101 37.1 


+68 09 53 


03-Mar-1996 


12:41:22 


3874 


V645 Cyg 


SOI 


26301850 


21 39 58.2 


+50 14 22 


06-Aug-1996 


07:29:22 


1912 


BD+46°3471 


C04 


54101787 


2152 34.1 


+47 13 44 


lO-May-1997 


17:20:17 


2804 


SVCep 


SOI 


28800703 


22 2133.0 


+73 40 24 


30-Aug-1996 


18:44:21 


3454 




P40 


56201203 


22 2133.1 


+73 40 27 


31 -May- 1997 


09:49:36 


364 


MWC 1080 


SOI 


26301659 


23 17 25.8 


+60 40 43 


06-Aug-1996 


04:43:40 


1912 




SOI 


28301459 


23 17 25.8 


+60 40 43 


25-Aug-1996 


23:15:57 


1912 



ISO-SWS was the Short Wavelength Spectrometer aboard 
ISO (2.38^5.2 /im). This instrument consisted of 4 sets of 12 
detectors. Each set of detectors covered a different wavelength 
region ([2.38:4.08], [4.08:12.0], [12.0:29.0] and [29.0:45.2] 
/im respectively). Not only the material of which the detectors 
were made differed from band to band, also the aperture sizes 
were larger at longer wavelengths (from 14"x 20" around 3 
yum to 20"x 33" around 30 fim). The data of different detec- 
tors and wavelength regions are calibrated independently. 

The source's spectrum was scanned in time, in a way that 
all 12 detectors of one set scanned through the same wave- 
length region twice. We only used the SWS data in which the 
whole spectrum was scanned (scanning mode AOTOl). In the 
standard pipeline reduction, dark currents (residual signal, not 
related to the actual astrophysical target) are subtracted from 
the raw data. The detector signal is corrected by multiplication 
with a responsitivity curve. Furthermore, a wavelength and flux 
calibration are carried out. We used the resulting product of 
the latest version of the off-line processing pipeline reduction 
(OLP 10). After this basic reduction, we manually removed 
bad data. This includes the suppression of deviating detectors 



(detectors which measure spectra that are too dissimilar with 
respect to the spectra observed by the other detectors operat- 
ing in the same wavelength region) and artefacts like glitches 
(spikes in intensity, spread out over several data points due to 
the detector's memory effects). The latter are caused by in- 
falling highly energetic interplanetary particles. A glitch that 
exceeds the noise level is easily recognizable, since it is only 
present in one of the two scan directions of the detector and its 
shape is characteristic. 

To diminish the scatter in the data, the calibrated flux lev- 
els in a wavelength interval are averaged and all detectors are 
scaled to this value. The noise distribution is strongly non- 
Gaussian and asymmetric, among other reasons because of 
glitches. Since the data will be rebinned and averaged, removal 
of the outliers is needed. Outliers with a deviation of more than 
three times the noise in the data were expelled. In the next step 
the data are rebinned at a spectral resolution of 150 for low flux 
sources (< 10 Jy at 10 yum) and 500 for high flux sources (> 10 
Jy at lOyum). 

Sometimes small differences in flux level are present be- 
tween the different wavelength regions in the spectrum. The fi- 
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nal phase of the data reduction is to correct for these jumps. 
Since most sources in the sample have low fluxes, jumps 
are primarily due to errors in the dark current subtraction. 
Therefore off'sets were applied to correct for the flux differ- 
ences between the wavelength intervals. The shifting was done 
by adding a constant term — which was computed to minimize 
the flux difference in the overlap region between two adjacent 
wavelength regions — to the flux. We took the interval between 
3.02 to 3.52 micron as the flux calibration reference and scaled 
all parts of the spectrum between 2.38 and 12.0 micron to that 
level. For the region [12.0:29.0] micron, the interval cover- 
ing [16.5:19.5] micron was used. We discarded the wavelength 
range longer than 29 micron, since the flux levels in this part of 
the spectrum were too low to be reliable for many of the sample 
stars. 

The wavelength region around 12 //m is difficult to inter- 
pret, since both the end of the previous and the beginning of 
the next wavelength range are very noisy and therefore un- 
reliable. For a few spectra, we found it useful to extend the 
wavelength region of the first interval with respect to the stan- 
dard pipeline reduction: from 12 to 12.5 jim. This was done 
for V376 Cas, HD 31648, Z CMa, HD 95881, HD 104237, HD 
139614, HD 142666, HD 150193, T CrA, WW Vul, LkHa 224, 
HD 200775 and SV Cep. Since we chose to shift the flux lev- 
els in [2.38:12.0]/im, and [12.0:29.0]//m separately, the overlap 
region between these two intervals is independent of this post- 
reduction flux level correction. In this way, we avoid changing 
the shape of spectrum around this wavelength area too much 
compared to the original — unshifted — data. Not combining 
the two intervals in a direct way also prevents the entry of sys- 
tematic errors in the shape and line flux of the broad, amor- 
phous 10 micron siUcate feature. We note that each of the nar- 
rower features at 3.3, 3.4, 3.5, 6.2, 7.7, 8.6, 11.2 and 11.3 mi- 
cron are included within one continuous SWS wavelength in- 
terval. Therefore our measurements of line flux and full width 
at half maximum of those features will be independent of the 
procedure followed for the alignment of the different wave- 
length ranges. 

For HD 97048, HD 100546, HD 104237, HD 142666, 
CD-42°1 1721, R CrA, T CrA and MWC 1080, two ISO-SWS 
spectra were available. We combined the spectra after the re- 
duction to increase the signal-to-noise ratio. For HD 97048 we 
decided to keep only one of the SWS spectra (OSN 61801318), 
because the second spectrum was affected by severe instrumen- 
tal artefacts. 

The ISO-PHT instrument was a photo-polarimeter, which 
had a spectroscopic mode (PHT-S). The spatial resolution 
of the instrument was 24"x24". The spectra consist of 127 
photometric points, separated into two wavelength intervals 
([2.5:4.9]/zm and [5.8:11.6]yL(m). For faint sources, PHT-S was 
able to obtain spectra with a higher S/N ratio than SWS, al- 
beit at much lower spectral resolution (~90). We retrieved the 
resulting spectra from the latest version of the off-line process- 
ing data reduction (OLP 10). For HD 97048, HD 104237, HD 
135344 and WW Vul, two or three PHT-S spectra were avail- 
able; we merged the spectra. 

The ISO-CAM instrument also had a spectrophotometric 
mode (CAM04). This mode made it possible to take a series of 



images at variable central wavelengths (between 5 and 15 fim). 
We extracted the spectrum by applying synthetic photometry 
to the images with a circular aperture with a diameter of 15" 
(comparable with the spatial resolution of SWS and CAM). 
Only 5 CAM04 spectra are included in this analysis. Again we 
used the results of the latest version of the off-line processing 
data reduction pipeline OLP 10. 

The resulting spectra of the data reduction are compiled 
in Fig.m Fig. Eland Fig.EJfor ISO-SWS, ISO-PHT and ISO- 
CAM, respectively. 

2.3. Measurement of the spectral features 

In our analysis, we focused on the infrared emission bands 
at 3.3, 6.2, 7.7, 8.6 and 11.2 micron, attributed to polycyclic 
aromatic hydrocarbons and the features at 3.4 and 3.5 mi- 
cron, linked to nanodiamonds. The amorphous silicate feature 
around 10 micron and the crystalline silicate feature at 1 1.3 mi- 
cron were also included in the analysis. The 1 1.2 and 1 1.3 mi- 
cron bands, when present, are blended. Due to the low spectral 
resolution, we cannot distinguish between the two. Therefore 
we consider the blend as one complex, which we call the 11 
micron feature. 

To characterize these features in the spectra, we measured 
the line flux, full width at half maximum (FWHM), peak wave- 
length, peak flux, continuum flux at the peak wavelength and 
the equivalent width (EW). 

We wrote an IDL procedure, measfeat.pro^, to manipu- 
late the spectra in order to determine the desired measurements. 
The IDL code allows us to indicate continuum points by hand. 
Based on these points a spline function is fitted to estimate 
the continuum. Integrating the continuum subtracted flux, the 
line flux and EW of the features were computed. We repeated 
the procedure, indicating an 'extreme' continuum, to estimate 
the systematic error introduced by the continuum determina- 
tion. It appears that the uncertainty in the Une flux varies from 
about ~10% for the PAH features and the 10 micron feature, 
up to ~30% for the 1 1 micron feature. Furthermore, the IDL 
program determined the peak position and peak flux by fitting 
a Gaussian function to the continuum-subtracted line profile. 
The peak flux was set equal to the maximum of the fit, the 
peak wavelength is the wavelength of this maximum. When 
the feature was in absorption, the peak flux was defined as the 
minimum of the fit. In the latter case the peak wavelength is 
the wavelength at which this minimum occurred. The value of 
the continuum flux at the peak wavelength is also recorded. 
Note that, because of the possible skewness of the features, the 
peak wavelength in general does not coincide with the centroid 
wavelength. This may lead to systematic differences between 
the values used in this paper and values in the literature. 

The FWHM was deduced by computing the width at half 
the peak flux of the smoothed line profile. The smoothing of 
the spectrum was done to diminish the influence of the noise in 
the spectrum on this measurement. 



^ this IDL procedure can be downloaded at the following 
website: http://www.ster.kuleuven.ac.be/~bram/ISO/measfeat.pro and 
_http://www.ster.kuleuven.ac.be/~bram/ISO/measfeat.README_ 
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Fig. 1. The reduced 2.4-15 /vm SWS spectra of the sample stars. 
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When a feature was not detected, we deduced upper lim- 
its for its line flux and peak flux, and we computed the con- 
tinuum flux at the expected central wavelength of the feature. 
The upper limits were computed in the following manner; we 
used the averaged specti'um of CD-42°11721 and HD 97048 
(both sources with fairly typical PAH spectra; Sect. 13. U as a 
template spectrum for the PAH emission, the spectrum of HD 
97048 for the 3.4 and 3.5 micron NAN emission and the spec- 
trum of HD 150193 as a template for the 10 micron amorphous 



silicate emission. For non-detections, we assumed a peak emis- 
sion of 5 times the noise on the spectrum at the theoretical cen- 
tral wavelength of the feature as an upper limit for the peak 
flux. Scaling the template spectra to this level, we extracted the 
upper limits for the line flux. 

This procedure was appUed to the SWS, PHT-S and CAM 
spectra. Since the spectral resolution of the PHT-S and CAM 
data is not even 100, only the line flux, EW, peak flux and con- 
tinuum flux at the theoretical central wavelength were deter- 
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Fig. 2. The reduced 2-5-1 1 .6 fj.m PHT-S spectra of the sample stars. 
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Fig. 3. The reduced 5-15 jum CAM spectra of the sample stars. 



mined for the latter; the other measurements are not meaningful 
in this case. 



2.4. Comparison of SWS, PHT-S and ground-based 
photometry 

For 24 sources, SWS as well as PHT-S spectra were available. 
We independently measured the desired quantities of the fea- 
tures, when present in both spectra. Comparing these two sets 
of values, we noticed that there is an offset between the two; the 
computed PHT-S line fluxes are systematically lower than their 



SWS counterparts. We attribute this difference to the fact that at 
2.8 times the Nyquist sample spacing, the PHT-S data are insuf- 
ficiently sampled in wavelength to recover the full spectral re- 
sponse. The result of this is that we will lose flux around sharp 
gradients in the spectra. As an example we show in Fig.l^the 
PAH 6.2 feature of HD 100546, as it appears in the SWS and 
PHT-S spectra. The line flux in the PHT-S spectrum is clearly 
underestimated. In Table|2lwe summarize the ratios of the line 
flux measured in the PHT-S spectrum over the line flux mea- 
sured in the SWS spectrum (PHT/SWS) for all detected fea- 
tures. Note that most of the ratios (81%) are smaller than unity. 
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Fig. 4. The PAH feature at 6.2 micron of HD 100546 in the 
SWS spectrum (dashed Une) and in the PHT-S spectrum (full 
line). 



which is consistent with our interpretation. The given errors are 
the statistical errors on the measurements. Ratios that are larger 
than one can be due to spatially extended emission or intrinsi- 
cally varying PAH emission. 

We calculated the mean ratio of the detected SWS line flux 
over the PHT-S line flux for each feature. The values of these 
quantities are 1.18, 1.34, 1.17, 1.16, 1.18, 1.08 and 0.58 for 
PAH 3.3, NAN 3.4, NAN 3.5, PAH 6.2, PAH 7.7, PAH 8.6 and 
the 1 1 micron feature respectively. The PHT-S measurements 
were adjusted by multiplying with these mean values. The fac- 
tors were also applied to correct the upper limits. 

Note that the amorphous 10 micron silicate feature is 
only partially sampled in the PHT-S spectra, since the wave- 
length coverage ends at 11.6 micron. Therefore it is difficult 
to determine the underlying continuum in the PHT-S spectra. 
Nevertheless, we indicated the continuum in a consistent man- 
ner. Similarly to the PAH features, we computed the line flux 
ratio of detected 10 micron emission features in both the SWS 
and the PHT-S spectrum, and computed the mean value. By 
applying this correction factor (0.86), we compensate for this 
systematical difference between the measurements of the SWS 
and PHT-S 10 micron feature. 

The comparison between the CAM and SWS spectra of HD 
97048 showed that the results are equal within a 15% range. No 
corrections for the CAM results were applied. We decided not 
to use the CAM spectra of HR 5999, HD 179218, HD 200775 
and HD 97048 in the analysis; the first two spectra are inhomo- 
geneously covered in wavelength, which leads to insufficiently 
sampled features. The CAM spectrum of HD 200775 is quite 
different from the SWS spectrum; the first displays PAH emis- 
sion, while the latter does not. We relate this to extended emis- 
sion, visible in the CAM spectrum because the (synthetic) aper- 
ture is larger than the aperture for the SWS spectrum at these 
wavelengths. The CAM spectrum of BD-h46°3471 is the only 
spectrum of this source at our disposal. We added the measure- 
ments to our data set, even though we noticed unusual features 



in the spectrum, which could be artefacts (e.g. the "absorption 
feature" at 10 micron). 

The ISO spectra were compared to space- and ground- 
based photometric measurements from the literature; we plot- 
ted the photometric points from the MSX catalog over the re- 
duced spectra. In general, the mean deviation between the spec- 
tra and the photometry is small and non-systematic. Sources 
where the anomaly exceeds the 20% level are HD 104237, 
HD 142666, MWC 297, BD-i-40°4124, V645 Cyg and MWC 
1080. Except for BD-i-40°4124, the photometric fluxes in these 
spectra are higher than the SWS or PHT-S fluxes. We attribute 
the differences to spatially extended emission and different 
beam/aperture sizes, although we cannot exclude the possibil- 
ity of intrinsic infrared variability of (some of) our sources. 

2.5. The final numbers 

After the measurements of the features, we condensed the data 
set into a final set of numbers: one measurement for each 
star and examined feature. The presence of the features in the 
HAEBE spectra is summarized in Table |3] The sources are 
sorted according to their IR-SED classification (see Sect. 12. 6> . 

When only a SWS, PHT-S or CAM spectrum was availaible 
for a source, those data were taken. For sources with a SWS 
and a PHT-S spectrum, three possibilities arise; the feature is 
detected in both spectra, the feature is undetected in both spec- 
tra or the feature is detected in only one of the two spectra. 
In the first case, the data are merged using a weighted aver- 
age, based on the statistical errors on the measurements. In the 
second case, the most stringent upper limit for the feature is 
taken. In the third case, when the measured line flux of the de- 
tected feature is lower than the upper limit derived from the 
other spectrum, the data are consistent and the values of the 
measurement are selected. In the inconsistent case (the upper 
limit is lower than the measured line flux), we set a conserva- 
tive upper limit for the feature. The latter only occurred in 4 
cases (out of 192 measured features). This behaviour might be 
either a result of the different aperture sizes used by SWS and 
PHT-S, in which case this could indicate the existence of spa- 
tially extended emission, or a first indication that the strength 
of PAH emission may vary with time in some of our sources. 
We indicated these measurements in Table |3l with V and V , 
according to detection in the PHT-S or SWS spectrum respec- 
tively. 

The computed line fluxes used in the analysis are listed in 
Table|4] Upper limits are indicated as well. Tablelsjsummarizes 
the continuum fluxes at the peak wavelength of each feature. 

2.6. The SED 

To characterize the spectral energy distribution (SED) of the 
sample sources, several quantities were determined, based on 
UV to millimetre (mm) photometry from the literature. The 
photometry consists of ANS and lUE ultraviolet data, ground- 
based photometry in the Walraven, Johnson/Cousins, and near- 
IR JHKLM photometric systems, IRAS and MSX data, and 
single-dish (sub)-mm photometry collected from the literature. 
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Table 2. The ratios of the line flux of the features measured in the PHT-S spectrum over the line flux measured in the SWS 
spectrum LFpht-s/LFsws- The given errors are the statistical errors on the measurements. 



Comparison PHT-S versus SWS 


Object 


PHT/SWS 3.3 


PHT/SWS 6.2 


PHT/SWS 7.7 


PHT/SWS 8.6 


PHT/SWS 11 


HD 34700 




<2.04 


<2.45 


<4.45 




HD 95881 




0.82 ± 0.06 


0.48 ± 0.03 


<1.33 


< 1.49 


HD 97048 


0.71 ±0.01 


1.06 ± 0.02 


0.97 ± 0.01 


0.70 ±0.01 


1.23 ±0.03 


HD 100453 


0.84 + 0.12 


0.59 ± 0.07 


1.00 ±0.05 


<0.63 




HD 100546 


0.88 + 0.04 


0.80 ± 0.01 


0.89 ±0.01 


0.87 ± 0.02 


2.55 ±0.10 


HD 135344 




<1.75 








HD 139614 




<2.83 


<2.21 


<1.52 


<6.67 


HD 141569 




<2.64 




<4.76 




HD 142527 


>0.87 


<1.14 


<1.41 


<1.91 


4.28 ± 0.28 


HD 142666 


>0.30 


3.11 ±0.24 


<2.00 




0.69 ±0.10 


HD 144432 




<1.23 






<0.97 


HR 5999 










1.55 ±0.19 


HD 150193 










< 1.47 


CD-42° 11721 


0.65 ± 0.01 


0.90 ± 0.01 


0.92 ± 0.01 


0.78 ± 0.01 


0.75 ± 0.01 


HD 163296 




0.91 ±0.08 






1.54 ±0.10 


HD 169142 


>1.61 


0.61 ± 0.05 


0.95 ± 0.03 


1.17 ±0.07 


1.92 ±0.23 


BD+40°4124 


>1.66 


0.51 ±0.01 


0.59 ± 0.01 


0.51 ±0.01 


1.00 ±0.03 



The effective temperature Teff of each source was determined 
from its spectral type listed in Table |6l using the calibrations 
from ISch midt-Kaler ( 1982). The visual extinction Ay and ab- 
solute luminosity L of the central star were computed by fitting 
and integrating a .Kurucz. C1991.) model for the stellar photo- 
sphere to the de-reddened photometry, adopting the distance 
estimates listed in Table|6l The ionising (>13.6 eV) luminosity 
Lion and the UV (2-13.6 eV) luminosity Luv of the source were 
computed using this Kurucz model fit. We also computed the 
IR-excess luminosity Lexc, using a spline fit to the infrared data, 
and the absorbed luminosity Labs of the source, which is the 
difference between the theoretical Kurucz model and the red- 
dened model. The IR-excess flux in the K band (2.2 ;um), at 60, 
850, and 1300 micron was deduced by subtracting the Kurucz 
model at these wavelengths from the observed fluxes. When no 
1300 micron photometry was available, we included 1 100 mi- 
cron measurements, which we then converted into 1300 micron 
flux values by multiplying with the average 1300/1 100 ratio de- 
duced from stars for which we had both measurements. The last 
parameter used to characterize the SED is the observed bolo- 
metric luminosity Lboi, which is the total integrated luminosity 
of the SED. The latter quantity is not corrected for extinction, 
as opposed to the stellar luminosity L. The stellar parameters 
used in this study are listed in Table|6l 

The sample sources were classified into different groups, 
based on the shape of their SED. The quantities used to char- 
acterize the infrared spectral energy distribution of HAEBEs 
are the ratio of Lnir (the integrated luminosity derived from 
the broad-band J, H, K, L and M photometry) and Lir (the 
corresponding quantity derived from IRAS 12, 25 and 60 mi- 
cron points), and t he non-color-c orrected IRAS ri21-r601 color 
Jvan Boekel et alJ i2003: Dullemond et alJ l2003t Ucke et alJ 
l2004l) . These parameters compare th e near-IR SED, which is 
expected to be similar in HAEBEs (Hillenbrand et al. 1992; 
[de Winter et al. 1995: Nattaet al. 2001) . to the mid-IR SED, 
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Fig. 5. Diagram based on Ivan Boekel et alJ ( 120031) . The ratio 
Lnir/Lir is plotted versus the non-color-corrected IRAS [12]- 
[60] color The dashed line represents Lnir/Lir - ([12]-[60])-i- 
0.9. Objects on the left-hand side of this empirical line are 
defined to be group II members (open diamonds), sources on 
the right side are either group I (filled diamonds) or group III 
(squares) members. 



where the major differences in SED shape occur. MOl Group I 
sources are stronger mid-IR emitters than group II sources. The 
luminosity ratio Lnir/Lir represents the strength of the near- 
IR compared to the mid-IR excess, which is lower for group 
I than group II sources. The shape of the mid-IR SED of a 
group I source is "double-peaked" compared to the SED of a 
group II member. The IRAS [12] -[60] color index provides a 
quantitative measure for this difference in SED shape. Group 
I sources are redder than their group II counterparts. We use 
the diagram given in Fig.lsjto distinguish between group I and 
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Table 3. The presence of the PAH features at 3.3, 6.2, 7.7 and 8.6 micron, the NAN features at 3.4 and 3.5 micron, the 1 1 micron 
complex (COMP 11) and the amorphous 10 micron silicate feature (Si 9.7) in the spectra of the sample stars. The sources are 
classified in 3 groups. ^J^. detection; - : no detection; n.s.: no spectrum available at these wavelengths. -y/° : tentative detection; 
y : tentative detection in the PHT-S spectrum; V' '■ feature undetected in the S WS spectrum; V : feature undetected in the PHT-S 
spectrum (see Sect. 12. 5> . For the amorphous silicate feature, E stands for emission and A for absorption. BD+40°4124, R CrA 
and LkHa 224 have not been classified based on the diagram in Fig.|5l which is indicated by the question mark. 



Group I/II vs IR emission bands. 


Object 


Group 


PAH 3.3 


NAN 3.4 


NAN 3.5 


PAH 6.2 


PAH 7.7 


PAH 8.6 


COMP 11 


Si 9.7 


V376 Cas 


1 


- 


- 


- 


- 


- 


- 


- 


E 


Ellas 3-1 


1 


V 


V 


V 


V 


V 


V 


V 


E 


AB Aur 


1 


- 


- 


- 


V 


V 


V 


V 


E 


HD 34282 


1 


- 


- 


- 


- 


V 


- 


- 


- 


HD 34700 


1 


- 


- 


- 


V 


V 


V 


- 


- 


HD 97048 


1 


V 


V 


V 


V 


V 


V 


V 


- 


HD 100453 


1 


V 


- 


- 


V 


V 


/ 


V 


- 


HD 100546 


1 


V 


/ 


- 


V 


V 


V 


V 


E 


HD 135344 


1 


- 


- 


- 


/ 


- 


- 


- 


- 


HD 139614 


1 


- 


- 


- 


V 


V 


V 


V 


- 


HD 142527 


1 


V 


- 


- 


- 


- 


- 


V 


E 


CD-42° 11721 


1 


V 


- 


- 


V 


V 


V 


V 


- 


HD 169142 


1 


v 


- 


- 


V 


V 


V 


V 


- 


TCrA 


1 


- 


- 


- 


- 


- 


- 


- 


E 


HD 179218 


1 


V 


- 


- 


V 


V 


V 


V 


E 


BD+40°4124 


1? 


/ 


V 


V 


V 


V 


V 


V 


- 


HD 200775 


1 


- 


- 


- 


V 


- 


- 


- 


E 


MWC 1080 


1 


V 


- 


- 


V 


V 


V 


V 


- 


VX Cas 


11 


V 


- 


- 


V 


V 


V 


- 


E 


HD 31648 


11 


- 


- 


- 


V 


- 


- 


- 


E 


UXOri 


11 


- 


- 


- 


- 


- 


- 


- 


E 


HD 35187 


11 


- 


- 


- 


- 


- 


- 


- 


- 


BFOri 


11 


- 


- 


- 


- 


- 


- 


- 


E 


RRTau 


11 


- 


- 


- 


V 


V 


V 


V 


- 


HD 95881 


11 


- 


- 


~ 


V 


V 


V 


V 


- 


HD 104237 


11 


- 


- 


- 


- 


- 


- 


- 


E 


HD 141569 


11 


- 


- 


- 


V 


V 


V 


- 


- 


HD 142666 


11 


V 


- 


- 


V 


- 


- 


V 


E 


HD 144432 


11 


- 


- 


- 


/ 


- 


- 


v 


E 


HR 5999 


11 


- 


- 


- 


- 


- 


- 


V 


E 


Wra 15-1484 


11 


V 


- 


- 


V 


V 


V 


V 


- 


HD 150193 


11 


- 


- 


- 


- 


- 


- 


V 


E 


AKSco 


11 


- 


- 


- 


- 


- 


- 


- 


E 


HD 163296 


11 


- 


- 


- 


V 


- 


- 


V 


E 


VVSer 


11 


- 


- 


- 


V" 


- 


- 


V 


E 


RCrA 


11? 


- 


- 


- 


- 


- 


- 


- 


E 


WWVul 


11 


- 


- 


- 


- 


- 


- 


- 


E 


LkHo- 224 


11? 


- 


- 


- 


- 


- 


- 


- 


- 


BD+46°3471 


11 


n.s. 


n.s. 


n.s. 


- 


- 


- 


- 


E 


SVCep 


11 


- 


- 


~ 


- 


- 


- 


- 


E 


ZCMa 


111 


- 


- 


~ 


- 


- 


- 


- 


A 


IRAS 12496-7650 


111 


- 


- 


- 


- 


- 


- 


- 


A 


MWC 297 


111 


V 


- 


V 


V 


V 


- 


V 


A 


LkHo- 225 


111 


- 


- 


- 


- 


- 


- 


- 


A 


PVCep 


111 


- 


- 


- 


- 


- 


- 


- 


A 


V645 Cyg 


111 


- 


- 


- 


- 


- 


- 


- 


A 



group II in the classification of MOl . The dashed line represents 
Lnir/^-ir - ([12] - [60]) + 0.9, which empirically provides the 
best separation between the two groups. 



Six of our sample stars display the amorphous 10 micron 
feature in absorption. These objects are believed to possess 
disks whose luminosity is dominated by viscous dissipation 
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Table 4. The line fluxes LF of the IR features included in this analysis. a(-b) represents a x 10 



-b 



LFoftheIR features. | 


Object 


Group 


PAH 3.3 


NAN 3.4 


NAN 3.5 


PAH 6.2 


PAH 7.7 


PAH 8.6 


COMP 1 1 


Si 9.7 






[W/m^] 


[W/m^] 


[W/m=] 


[W/m^] 


[W/m^] 


[W/m^] 


[W/m^] 


[W/m=] 


V376 Cas 


I 


<2.38(-15) 


<2.82(-15) 


<2.26(-I5) 


<1.62(-14) 


<2.34(-14) 


<9.64(-I5) 


<8.88(-15) 


5.39(-13) 


Elias 3-1 


I 


7.28(-15) 


1.13(-14) 


1.38(-14) 


4.49(-14) 


1.23(-I3) 


3.7I(-14) 


7.40(-I4) 


9.11(-I3) 


AB Aui- 


I 


<5.33(-15) 


<3.26(-15) 


<2.61(-I5) 


3.99(-I4) 


5.63(-14) 


3.68(-I4) 


4.61(-I4) 


I.47(-I2) 


HD 34282 


I 


<2.81(-15) 


<3.03(-I5) 


<2.43(-I5) 


<2.43(-14) 


2.74(-I4) 


<I.16(-I4) 


<I.08(-I4) 


<I.06(-I3) 


HD 34700 


I 


<2.81(-15) 


<3.45(-15) 


<2.49(-15) 


8.75(-15) 


I.46(-I4) 


3.28(-I5) 


<I.28(-I5) 


<9.74(-I5) 


HD 97048 


I 


1.85(-14) 


4.71(-14) 


8.03(-14) 


8.83(-14) 


2.40(-I3) 


4.52(-I4) 


3.97(-14) 


<8.95(-I5) 


HD 100453 


I 


8.99(-15) 


<3.86(-15) 


<3.I0(-I5) 


2.69(-I4) 


2.70(-I4) 


<I.65(-I4) 


7.83(-I5) 


<3.34(-I4) 


HD 100546 


I 


2.40(-I4) 


8.80(-15) 


<I.74(-I5) 


I.2I(-I3) 


3.22(-I3) 


5.6I(-I4) 


3.89(-I3) 


3.I8(-I2) 


HD 135344 


I 


<3.06(-15) 


<3.03(-15) 


<2.43(-I5) 


I.09(-14) 


<9.I6(-I5) 


<3.6I(-I5) 


<I.56(-I5) 


<2.4I(-I4) 


HD 139614 


I 


<2.83(-15) 


<3.22(-I5) 


<2.59(-I5) 


6.07(-I5) 


I.40(-I4) 


8.32(-I5) 


1.23(-I5) 


<2.9I(-I4) 


HD 142527 


I 


7.21(-15) 


<3.56(-15) 


<2.85(-I5) 


<1.69(-I4) 


<2.4I(-I4) 


<9.90(-I5) 


2.09(-I4) 


5.52(-I3) 


CD-42° 1 1721 


I 


6.87(-I4) 


<4.66(-I5) 


<3.73(-I5) 


4.89(-I3) 


I.34(-I2) 


4.13(-I3) 


3.54(-I3) 


<9.08(-I4) 


HD 169142 


I 


<7.38(-15) 


<2.96(-I5) 


<2.37(-I5) 


1.93(-I4) 


5.68(-I4) 


2.55(-I4) 


8.I8(-I5) 


<2.I9(-I4) 


TCrA 


I 


<1.43(-15) 


<1.72(-I5) 


<1.38(-15) 


<3.I0(-I5) 


<7.I7(-I5) 


<3.75(-I5) 


<3.24(-I5) 


I.92(-I3) 


HD 179218 


I 


1.52(-I4) 


<2.53(-15) 


<2.03(-I5) 


7.07(-I4) 


I.42(-I3) 


2.66(-I4) 


1.39(-I3) 


I.03(-I2) 


BD+40°4124 


I? 


<1.20(-14) 


4.10(-15) 


4.13(-15) 


5.38(-14) 


9.84(-I4) 


3.32(-I4) 


5.37(-I4) 


<4.70(-I4) 


HD 200775 


I 


<3.08(-15) 


<3.25(-15) 


<2.60(-15) 


3.35(-I4) 


<2.05(-I4) 


<9.I2(-I5) 


<I.I9(-I4) 


3.36(-I3) 


MWC 1080 


I 


I.I4(-14) 


<3.84(-I5) 


<3.08(-I5) 


3.64(-14) 


9.09(-I4) 


2.90(-I4) 


5.5I(-I4) 


<I.I6(-I3) 


VX Cas 
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3.00(-I5) 


<2.I4(-I5) 


<I.50(-I5) 


3.17(-15) 


6.25(-15) 


2.7I(-15) 


<7.28(-I6) 


3.9I(-I4) 


HD 31648 


II 


<2.4I(-I5) 


<3.25(-I5) 


<2.61(-I5) 


I.59(-14) 


<2.I9(-I4) 


<I.02(-I4) 


<9.59(-I5) 


5.85(-I3) 


UXOri 


II 


<1.4I(-15) 


<I.76(-15) 


<I.23(-I5) 


<6.53(-16) 


<2.32(-I5) 


<I.57(-I5) 


<7.20(-I6) 


I.32(-I3) 


HD 35187 


II 


<5.4I(-15) 


<3.35(-I5) 


<2.69(-I5) 


<1.4I(-14) 


<2.36(-I4) 


<I.06(-I4) 


<I.I5(-I4) 


<I.09(-I3) 


BF Oi-i 


II 


<1.39(-15) 


<1.74(-I5) 


<I.21(-I5) 


<7.38(-I6) 


<2.35(-I5) 


<I.26(-I5) 


<6.I0(-I6) 


5.II(-I4) 


RRTau 


II 


<1.62(-15) 


<2.01(-15) 


<I.41(-15) 


I.07(-I4) 


2.38(-I4) 


I.0I(-I4) 


2.53(-I5) 


<5.43(-I5) 


HD 95881 


II 


<3.58(-15) 


<2.32(-I5) 


<I.86(-I5) 


2.15(-14) 


I.96(-I4) 


8.12(-I5) 


4.05(-I5) 


<4.I2(-I4) 


HD 104237 


II 


<2.54(-15) 


<2.98(-15) 


<2.39(-I5) 


<2.00(-I5) 


<6.52(-I5) 


<2.67(-I5) 


<I.56(-I5) 


7.32(-I3) 


HD 141569 


II 


<2.42(-15) 


<3.01(-I5) 


<2.I1(-I5) 


6.91(-15) 


I.62(-I4) 


2.94(-I5) 


<8.02(-I6) 


<7.8I(-I5) 


HD 142666 


II 


2.17(-I5) 


<I.73(-I5) 


<I.38(-I5) 


2.37(-I4) 


<1.30(-14) 


<5.37(-I5) 


5.42(-I5) 


3.85(-I3) 


HD 144432 


II 


<1.75(-15) 


<1.94(-I5) 


<I.55(-15) 


8.98(-15) 


<I.I5(-I4) 


<4.75(-I5) 


3.54(-I5) 


3.88(-I3) 


HR 5999 


II 


<4.65(-15) 


<4.55(-I5) 


<3.65(-I5) 


<7.54(-I5) 


<1.74(-I4) 


<7.I7(-I5) 


9.95(-I5) 


2.78(-I3) 


Wra 15-1484 


II 


5.57(-15) 


<4.60(-15) 


<3.69(-I5) 


4.76(-14) 


I.I6(-I3) 


2.43(-I4) 


2.36(-I4) 


<8.36(-I4) 


HD 150193 


II 


<3.06(-15) 


<3.80(-I5) 


<2.66(-I5) 


<I.8I(-I5) 


<4.96(-I5) 


<3.64(-I5) 


8.29(-I5) 


I.12(-I2) 


AKSco 


II 


<2.90(-15) 


<3.60(-15) 


<2.52(-I5) 


<I.65(-15) 


<5.40(-I5) 


<2.22(-I5) 


<I.I3(-I5) 


7.88(-I4) 


HD 163296 


II 


<3.04(-15) 


<3.I5(-15) 


<2.53(-15) 


2.61(-14) 


<I.35(-I4) 


<5.57(-I5) 


1.49(-I4) 


8.41(-I3) 


VVSer 


II 


<2.63(-15) 


<3.27(-I5) 


<2.29(-I5) 


I.06(-14) 


<3.30(-I5) 


<I.67(-I5) 


2.04(-I5) 


8.83(-I4) 


RCrA 


II? 


<8.35(-15) 


<9.63(-15) 


<7.72(-I5) 


<I.85(-14) 


<4.58(-I4) 


<2.25(-I4) 


<I.28(-I4) 


I.67(-I2) 


WW Vul 


II 


<2.00(-15) 


<2.49(-I5) 


<I.74(-15) 


<I.43(-I5) 


<3.69(-I5) 


<I.5I(-I5) 


<I.09(-I5) 


I.27(-I3) 


LkHo- 224 


II? 


<3.14(-15) 


<3.59(-15) 


<2.88(-I5) 


<1.42(-14) 


<2.4I(-I4) 


<I.I3(-I4) 


<I.0I(-I4) 


<9.30(-I4) 


BD-i-46°3471 


II 








<5.90(-15) 


<1.42(-I4) 


<7.28(-I5) 


<4.75(-I5) 


3.I9(-I4) 


SVCep 


II 


<1.67(-15) 


<2.I9(-15) 


<I.76(-I5) 


<7.58(-16) 


<2.93(-I5) 


<I.20(-I5) 


<9.75(-I6) 


2.35(-I3) 


ZCMa 


III 


<6.98(-15) 


<8.29(-I5) 


<6.65(-15) 


<4.86(-14) 


<1.22(-13) 


<5.48(-I4) 


<4.60(-I4) 


-3.50(-I2) 


IRAS 12496-7650 


III 


<1.09(-14) 


<1.35(-I4) 


<9.50(-15) 


<5.41(-I5) 


<1.51(-14) 


<6.21(-15) 


<2.75(-I5) 


-4.98(-I3) 


MWC 297 


III 


1.93(-I4) 


<1.06(-I4) 


2.73(-14) 


I.54(-13) 


I.89(-I3) 


<2.98(-I4) 


4.88(-I4) 


-9.47(-I3) 


LkHo- 225 


III 


<2.25(-15) 


<2.64(-15) 


<2.12(-I5) 


<I.I5(-I4) 


<I.77(-I4) 


<7.97(-I5) 


<9.30(-I5) 


-I.80(-12) 


PVCep 


III 


<2.67(-15) 


<3.05(-I5) 


<2.45(-I5) 


<2.76(-15) 


<1.10(-14) 


<4.52(-I5) 


<2.24(-I5) 


-4.99(-I3) 


V645 Cyg 


III 


<3.65(-15) 


<4.54(-I5) 


<3.64(-I5) 


<2.00(-14) 


<3.64(-14) 


<1.63(-14) 


<2.04(-I4) 


-2.33(-I2) 



of energy due to accretion, are deeply embedded systems and 
hence are fundamentally diff'erent from the other sample stars. 
We therefore classify them in a difi'erent group: group III. 

BD+40°4124, R CrA and LkRa 224 have not been classi- 
fied based on their appearance in Fig.|5] Confusion with back- 
ground sources in the IRAS photometry prohibited us from 
deriving the quantities needed to plot these objects in the di- 
agram. BD+40°4124 has been classified as a group I source, 
because its SED resembles the SED of HD 200775. R CrA and 
LkHa 224 are both UX Orioni s stars according to the defini- 
tion of IPuUemond et alJ 1120031) : both sources have a spectral 
type later than B9 and optical variations larger than 1 mag on 
timescales of days to weeks. Therefore we classified them as 
group II members. 

The SEDs of all sample sources are displayed in Fig. |6l 
The group to which the objects belong is indicated as well. The 
SEDs are overplotted with the ISO spectra. Squares with error 
bars indicate photometric measurements from the literature, ar- 
rows upper limits. 



3. Analysis of the ISO spectra 

3.1. The PAH features 

First we investigated the general correlations of the different 
PAH features in our spectrum. We compared the line fluxes 
(LF) of the different features, which we normalised by the con- 
tinuum flux at the peak wavelength (CF), to correct for the 
distance-to-the-star dependence of the line flux. Fig. shows 
the ratio LF/CF for the RAH 7.7 micron feature versus PAH 
6.2 micron; Fig.|Slis a similar plot for the ratio LF/CF of PAH 
6.2 micron versus PAH 3.3 micron. The correlation coeflicient 
T] in the first plot is 0.90, while the data in the second plot are 
less correlated [j] = 0.58). The other PAH features display com- 
parable behaviour. It is striking though that PAH features that 
are linked to the CC bonds (6.2, 7.7 micron) correlate better 
with each other, as well as the features linked to the CH bonds 
(3.3, 7.7, 8.6 micron). However, as a first order approximation 
we can state that the LF/CF ratios correlate fairly well over a 
range of a few orders of magnitude. 



12 



Acke & v.d. Ancker: Disks Around HAEBE Stars 



Table 5. The continuum flux at the peak wavelength of the IR features included in this analysis. When the amorphous 10 micron 
silicate band was present, the PAH features are superimposed on the band. In this case (flagged by *), the silicate feature is 
regarded as being part of the underlying continuum of the PAH features. Hence the indicated continuum flux of the PAH band is 
the flux of the silicate feature at the central wavelength of that PAH band. 



CF at the peak wavelength of the IR features. 


Object 


Group 


PAH 3.3 


NAN 3.4 


NAN 3.5 


PAH 6.2 


PAH 7.7 


PAH 8.6 


COMPll 


Si 9.7 






[Jy] 


[ly] 


[Jy] 


[Jy] 


[Jy] 


[Jy] 


[Jy] 


[Jy] 


V376 Cas 


I 


5.5 


6.0 


6.3 


10.0 


11.3* 


17.1* 


23.2* 


13.7 


Ellas 3-1 


I 


3.1 


3.3 


3.5 


4.0 


8.0* 


12.4* 


32.7* 


19.1 


AB Aur 


I 


13.1 


13.1 


13.0 


9.1 


10.0* 


15.1* 


24.8* 


12.8 


HD 34282 


I 


0.7 


0.8 


0.7 


1.0 


0.9 


0.8 


0.8 


0.7 


HD 34700 


I 


0.5 


0.5 


0.3 


0.4 


0.7 


0.7 


0.8 


0.7 


HD 97048 


I 


2.0 


2.0 


2.0 


1.1 


3.2 


4.0 


6.4 


4.7 


HD 100453 


I 


4.8 


4.9 


4.9 


3.7 


4.4 


4.2 


6.3 


4.7 


HD 100546 


I 


6.0 


5.7 


5.7 


7.0 


17.3* 


28.8* 


62.5* 


28.4 


HD 135344 


I 


3.4 


3.4 


3.4 


2.1 


1.7 


1.4 


1.0 


1.1 


HD 139614 


I 


1.4 


1.3 


1.3 


1.0 


1.4 


1.5 


2.6 


2.4 


HD 142527 


I 


6.4 


6.4 


6.4 


7.0 


1.1' 


10.7* 


11.9* 


7.1 


CD-42ni721 


I 


23.6 


25.0 


25.7 


41.4 


44.3 


43.7 


45.9 


39.5 


HD 169142 


I 


1.8 


1.7 


1.7 


1.7 


1.6 


1.6 


2.0 


1.6 


TCrA 


I 


2.2 


2.2 


2.3 


3.0 


3.3* 


4.4* 


8.3* 


5.9 


HD 179218 


I 


3.4 


3.4 


3.5 


5.5 


9.3* 


14.2* 


18.3* 


10.5 


BD+40°4124 


I? 


5.5 


5.6 


5.7 


7.5 


9.0 


9.6 


8.9 


8.1 


HD 200775 


I 


11.8 


11.7 


11.8 


9.4 


9.8- 


11.4* 


9.3* 


8.2 


MWC 1080 


I 


11.2 


11.6 


11.7 


14.5 


14.9 


14.6 


16.2 


13.6 


VXCas 


II 


0.4 


0.4 


0.4 


0.4 


0.5* 


0.8* 


1.2* 


0.7 


HD 31648 


II 


4.0 


4.0 


4.0 


4.6 


6.1* 


8.7* 


11.3* 


6.1 


UXOri 


II 


0.9 


0.9 


0.9 


0.8 


0.9* 


2.0* 


2.7* 


1.5 


HD 35187 


II 


2.4 


2.4 


2.4 


1.9 


2.6 


3.8 


5.4 


4.9 


BFOri 


II 


0.4 


0.4 


0.4 


0.5 


0.6* 


1.1* 


1.4* 


0.8 


RRTau 


II 


0.9 


0.9 


0.9 


0.9 


1.1 


1.2 


1.6 


1.3 


HD 95881 


II 


4.0 


3.9 


3.9 


4.5 


5.5 


6.9 


7.8 


8.0 


HD 104237 


II 


10.8 


10.8 


10.8 


10.8 


10.6* 


15.9* 


22.4* 


9.6 


HD 141569 


II 


0.7 


0.7 


0.6 


0.4 


0.5 


0.5 


0.5 


0.5 


HD 142666 


II 


2.7 


2.7 


2.7 


0.6 


3.7* 


5.5* 


6.3* 


2.6 


HD 144432 


II 


2.8 


2.8 


2.7 


2.1 


2.3* 


4.8* 


7.4* 


4.1 


HR 5999 


II 


17.9 


18.0 


18.1 


13.6 


13.9* 


15.0* 


16.9* 


14.7 


Wra 15-1484 


II 


6.6 


7.1 


7.4 


14.4 


17.3 


18.4 


21.0 


19.8 


HD 150193 


II 


4.8 


4.7 


4.7 


4.8 


5.5* 


13.3* 


18.7* 


8.8 


AKSco 


II 


1.5 


1.5 


1.5 


1.2 


1.2* 


1.9* 


2.8* 


2.0 


HD 163296 


II 


9.8 


9.7 


9.6 


7.2 


8.8* 


14.0* 


20.3* 


12.2 


VVSer 


II 


2.5 


2.5 


2.5 


2.6 


3.1* 


3.8* 


4.4* 


3.4 


RCrA 


II? 


54.2 


57.0 


59.2 


64.6 


58.3* 


79.8* 


126.3* 


86.1 


WWVul 


II 


0.9 


0.9 


0.8 


0.7 


0.8* 


1.5* 


2.1* 


1.3 


LkHo- 224 


II? 


3.9 


4.1 


4.2 


7.3 


9.0 


10.3 


12.4 


11.4 


BD+46°3471 


II 








0.6 


0.9* 


1.0* 


1.6* 


1.0 


SVCep 


II 


0.6 


0.7 


0.7 


0.8 


1.1* 


2.9* 


4.6* 


2.7 


ZCMa 


III 


40.1 


43.3 


45.2 


86.9 


101.7 


84.3 


101.9 


114.7 


IRAS 12496-7650 


III 


7.9 


9.1 


9.7 


17.7 


21.8 


21.8 


33.3 


30.2 


MWC 297 


III 


65.8 


68.3 


71.5 


104.6 


118.9 


123.7 


119.7 


123.3 


LkHo- 225 


III 


1.3 


1.5 


1.7 


12.6 


24.5 


23.8 


26.2 


36.8 


PVCep 


III 


1.9 


2.2 


2.4 


7.4 


9.6 


8.5 


10.6 


12.8 


V645 Cyg 


III 


8.3 


9.6 


10.7 


45.2 


65.7 


65.4 


87.5 


91.0 



Supposing that all HAEBEs have chemically identical PAH 
molecules, the PAH spectra of this group of stars should be uni- 
form. Two effects that could influence the shape of this 'unique' 
spectrum are the UV field of the central star and interstellar ex- 
tinction. If we presume that the wavelength dependence of the 



first does not vary much from source to source, the dominant 
source of deviations from a simple linear correlation between 
the PAH feature-strengths would be extinction. This will have a 
stronger effect on PAH features that are close to the amorphous 
10 micron silicate feature. In Fig.|5]the logarithm of the ratio of 



Table 6. The stellar parameters for the sample sources. The distance d which we used in the analysis and spectral type (Sp.T.) are given. The visual extinction Ay, the logarithm 
of the effective temperature log T^ff, of the stellar luminosity log L, of the stellar UV luminosity log Luv and of the observed bolometric luminosity log Lboi are presented as 
well. The latter quantities are estimated based on the photometric measurements. Furthermore, the ratio of the IR excess luminosity over the stellar luminosity is computed. 
Ten of the sample sources have a ratio larger than unity, hence are not consistent with a passive reprocessing disk. The l ast 4 c olumns contain the IR excess fluxes in K (2.2 
jum), at 60, 850 and 1300 micron. * These stars are UX Orionis objects according to the definition o f Dullemond et alJ ( 120031) . * Probable confusion with the nearby source 
CrA IRS7. 











The stellar parameters of the sample s 


stars. 














Object 


Region 


d 


Ref. 


Sp.T 


Ref. 


Av 


log Teff 


logL 


logL„v 


log Lboi 


Lexc/L 


AK 


A 60 


A 850 


A 1300 






[pc] 








[mag] 


log[K] 


log[Lo] 


log[Lo] 


log[Lo] 




[mag] 


[mag] 


[mag] 


[mag] 


group I 
































V376 Cas 


CasRl 


630 


S89 


B5e 


C79 


4.40 


4.188 


1.63 


1.60 


2.57 


9.35 


5.65 


16.45 


15.48 


14.61 


Elias3-1 


LI 495 


160 


K94 


A6:e 


Z94 


4.05 


3.904 


-0.14 


-0.33 


1.35 


35.69 


5.56 


14.89 


15.92 


15.78 


AB Aur 


L1519 


140 


V98 


AOVe+sh 


B93 


0.50 


3.979 


1.68 


1.56 


1.70 


0.51 


2.28 


11.33 


11.08 


10.67 


HD 34282* 




400 


P03 


A3Vne 


MOO 


0.28 


3.941 


1.27 


1.11 


1.37 


0.53 


2.05 


11.68 


14.04 


- 


HD 34700 


Orion OB la 


340 


D99 


GOIV+GOIVe 


A03 


0.00 


3.774 


1.31 


0.91 


1.42 


0.27 


0.08 


10.21 


- 


9.96 


HD 97048 


Cedlll 


180 


V98 


B9.5Ve+sh 


W87 


1.26 


4.000 


1.64 


1.53 


1.41 


0.35 


1.23 


11.52 


- 


13.04 


HD 100453 




112 


HIP 


A9Ve 


H75 


0.02 


3.869 


0.90 


0.66 


1.09 


0.57 


1.27 


10.71 


- 


- 


HD 100546 


Sco OB2-4? 


103 


V98 


B9Vne 


H75 


0.26 


4.021 


1.51 


1.42 


1.62 


0.58 


1.49 


11.77 


- 


- 


HD 135344 


Sco OB2-3 


140 


D99 


F4Ve 


D97 


0.31 


3.819 


0.91 


0.59 


1.00 


0.44 


1.23 


10.23 


12.04 


11.64 


HD 139614 


Sco OB2-3 


140 


D99 


A7Ve 


D97 


0.09 


3.895 


0.91 


0.70 


1.03 


0.40 


0.98 


10.57 


12.72 


12.80 


HD 142527 


Sco OB2-2 


145 


D99 


F7IIIe 


V98 


0.64 


3.796 


1.18 


0.82 


1.34 


0.82 


1.28 


11.08 


13.47 


13.18 


CD-42° 11721 




400 


D90 


BOIVep 


S90 


5.08 


4.470 


3.94 


3.92 


3.02 


0.13 


2.83 


14.87 


12.19 


12.07 


HD 169142 


Sco OB 2-1 


145 


D99 


A5Ve 


D97 


0.43 


3.914 


1.16 


0.98 


1.12 


0.26 


0.76 


10.66 


12.38 


12.18 


TCrA 


NGC 6729 


130 


M81 


FOe 


F84 


2.45 


3.857 


-0.18 


-0.44 


0.88 


12.21 


3.86 


14.38 


<13.29 


14.99 


HD 179218 


L693 


240 


V98 


B9e 


S66 


0.54 


4.021 


2.00 


1.90 


1.88 


0.28 


1.27 


10.53 


- 


10.93 


BD+40°4124 


CygRl 


980 


S91 


B2Ve 


H95 


3.01 


4.342 


3.77 


3.76 


2.49 


0.05 


2.28 


- 


12.72 


13.34 


HD 200775 


CepR2 


440 


W81 


B2.5IVe 


R95 


1.90 


4.291 


3.82 


3.80 


2.84 


0.07 


1.44 


12.49 


9.48 


8.17 


MWC 1080 


L1238 


2200 


L88 


BOe 


C79 


5.27 


4.477 


5.26 


5.24 


3.75 


0.03 


2.98 


12.41 


12.60 


12.44 


group II 
































VX Cas* 


CasRl 


630 


S89 


AlVe+sh 


G98 


0.81 


3.965 


1.50 


1.37 


1.44 


0.47 


1.61 


10.04 


- 


<11.00 


HD 31648 




130 


V98 


A3Ve 


J91 


0.25 


3.941 


1.14 


0.98 


1.22 


0.42 


1.79 


9.63 


12.76 


12.87 


UX Ori* 


Orion OB la 


340 


D99 


A4IVe 


MOO 


1.91 


3.925 


1.68 


1.51 


1.25 


0.16 


0.91 


8.73 


- 


10.49 


HD 35187 


L1559 


150 


D98 


A2Ve+A7Ve 


D98 


0.71 


3.953 


1.44 


1.29 


1.23 


0.15 


0.99 


8.91 


10.20 


9.99 


BF Ori* 


OriOBlc 


510 


D99 


A2IVev 


MOO 


0.87 


3.954 


1.53 


1.39 


1.38 


0.30 


1.56 


9.68 


- 


10.51 


RR Tau* 


(L1553) 


160 


K94 


A3-5e 


F84 


1.04 


3.927 


0.31 


0.14 


0.42 


0.96 


2.73 


11.06 


- 


<12.15 


HD 95881 


Sco OB2-4? 


118 


D99 


A2III/IVe 


H75 


0.25 


3.954 


0.84 


0.69 


0.98 


0.62 


2.15 


8.06 


- 


- 
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The stellar parameters of the sample stars (continued). 


Object 


Region 


d 


Ref. 


Sp.T 


Ref. 


Ak 


log Teff 


logL 


log L„v 


logLboi 


Lexc/L 


AK 


60 


850 


1300 






[pc] 










log[K] 


log[Lo] 


log[Lo] 


log[Lo] 




[mag] 


[mag] 


[mag] 


[mag] 


HD 104237 


Cha III 


116 


V98 


A4IVe+sh 


V98 


0.29 


3.925 


1.54 


1.36 


1.53 


0.23 


1.18 


8.54 


- 


10.02 


HD 141569 


(LI 69) 


99 


V98 


AOVe 


D97 


0.37 


3.979 


1.28 


1.16 


1.09 


0.01 


-0.07 


8.18 


11.16 


8.21 


HD 142666* 


Sco OB2-2 


145 


D99 


A8Ve 


D97 


0.93 


3.880 


1.13 


0.91 


1.03 


0.33 


1.21 


8.96 


11.54 


11.51 


HD 144432 


Sco OB2-2 


145 


D99 


A9IVev 


MOO 


0.17 


3.866 


1.01 


0.76 


1.13 


0.45 


1.40 


8.89 


10.41 


10.37 


HR 5999* 


Lupus 3 


210 


V98 


A5-7III/IVe+sh 


T89 


0.49 


3.899 


1.94 


1.74 


1.97 


0.46 


1.63 


7.98 


9.08 


<9.19 


Wra 15-1484 




750 


L89 


BO:[e] 


D98b 


3.07 


4.477 


2.99 


2.97 


2.45 


0.31 


4.45 


13.60 


- 


- 


HD 150193 


Sco OB2-2 


150 


V98 


Al-3Ve 


G98 


1.49 


3.953 


1.38 


1.23 


1.19 


0.45 


1.82 


8.98 


10.57 


10.43 


AK Sco* 




150 


V98 


F5+F5IVe 


A89 


0.62 


3.809 


0.95 


0.62 


0.88 


0.21 


0.64 


8.69 


10.07 


- 


HD 163296 




122 


V98 


A3Ve 


G98 


0.09 


3.941 


1.38 


1.22 


1.51 


0.46 


1.92 


9.98 


13.00 


12.97 


VV Ser* 


Serpens 


330 


D91 


A2IIIe 


V 


2.67 


3.954 


1.27 


1.13 


1.36 


1.37 


3.13 


10.08 


- 


<11.39 


R CiA* 


NGC 6729 


130 


M81 


Al-F7ev 


V98 


1.33 


3.857 


-0.19 


-0.45 


3.45 


4675.32* 


6.56 


<17.10 


<15.68 


<14.89 


WW Vul* 


VulRl 


440 


V81 


A4IVA'e+sh 


G98 


1.18 


3.979 


1.44 


1.31 


1.29 


0.50 


2.19 


9.71 


- 


11.28 


LkHa 224* 


CygRl 


980 


S91 


A7e 


V99 


2.98 


3.895 


2.06 


1.85 


2.31 


1.97 


3.15 


- 


< 14.06 


- 


BD+46°3471 


IC5146 


1200 


H02 


A0.5IIIe 


G98 


1.27 


3.991 


2.60 


2.49 


2.43 


0.50 


2.47 


9.44 


- 


<11.17 


SV Cep* 


CepR2 


440 


W81 


A0-2Ve 


G98 


0.68 


3.965 


1.15 


1.01 


1.39 


1.29 


2.49 


10.95 


- 


11.41 


group III 
































ZCMa 


CMaRI 


1050 


S99 


Be 


V98 


4.46 


4.477 


5.15 


5.12 


3.71 


0.05 


2.57 


11.88 


11.98 


11.85 


IRAS 12496-7650 


Cha II 


180 


W97 


Ae 


H91 


14.00 


3.927 


2.52 


2.35 


1.60 


0.14 


0.51 


9.16 


- 


10.71 


MWC 297 


(L5I5) 


250 


D97c 


B1.5Ve 


D97c 


7.73 


4.375 


4.01 


4.00 


2.69 


0.07 


2.61 


12.11 


10.31 


10.76 


LkHo- 225 


CygRl 


980 


S91 


A-Fe 


H95 


7.32 


3.857 


2.62 


2.36 


3.25 


4.64 


2.05 


13.41 


14.19 


13.79 


PVCep 


CepR2 


440 


W81 


A5:e 


C81 


7.01 


3.914 


0.94 


0.75 


1.91 


11.29 


4.33 


14.24 


15.98 


15.81 


V645 Cyg 


Cyg OB7 


3500 


G86 


07e 


V02 


4.15 


4.580 


4.54 


4.41 


4.60 


1.25 


6.23 


16.99 


17.53 


17.26 
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Fig. 6. The spectral energy distributions of the sample stars. The squares are the observed photometric data from the literature 
with error bars. The solid line is the reddened Kurucz model for the stellar photosphere, fitted to the measured UV-optical 
photometry. The ISO spectra are plotted in the SEDs as well (noisy solid line). 



the line fluxes of the PAH feature at 8.6 micron over the feature terstellar extinction law o nFluks et alJ ( Il994h to the line fluxes 

at 6.2 micron, log(LF^,(,/LF(,2), is plotted versus the logarithm of MWC 1080 according to the formula 

of the ratio of the line fluxes of the 3.3 micro n PAH feature ove r 

the 6.2 feature, logiLF^.^ILFe.!} (following lSonv et al.l20oil) . 

To include the extinction eff'ect into the plot, we appUed the in- LF'(x) - LF^^'^ '"**" x 10-^ ^ 
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Fig. 7. The LFICF ratio of PAH 7.7 micron versus PAH 6.2 
micron. The filled diamonds with error bars indicate detected 
features, arrows represent upper limits. 
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Fig. 8. Similar plot as Fig. IT] The LFjCF ratio of PAH 6.2 
micron versus PAH 3.3 micron. 

in which x represents the actual visual extinction (in mag) and 
A^jAy the interstellar extinction coefficient at the central wave- 
length A of feature j. The straight line in Fig.|9]represents the 
ratios of extinction corrected line fluxes for a wide range of vi- 
sual extinction. If interstellar extinction would be the only ef- 
fect, all measurements should lie within the observational error 
bars from this line. This is not the case; sources that deviate sig- 
nificantly are Elias 3-1, AB Aur, RR Tau and CD-42°11721. A 
similar plot based on the PAH 7.7 micron feature instead of the 
8.6 micron feature (not in the paper), leads to the same conclu- 
sion: the spread cannot be induced by extinction only. Hence, 
the PAH spectra of HAEBE stars are not identical, indicating 
differences between PAH molecules in different systems. This 
spectral variety of PAHs can have several causes. 

1 . As described in the introduction, the different PAH emis- 
sion bands studied in this analysis are linked to bending 
and stretching modes in the different bonds of the PAH 



CO 



CD 
CO 



O 




log(LF_3.3/LF_6.2) 

Fig. 9. The logarithm of the ratio of the line fluxes of the 
PAH features at 8.6 and 6.2 micron \Qg{LF%,^ j LF (,2) versus the 
logarithm of the ratio of the line fluxes of the PAH features 
at 3.3 and 6.2 micron \og{ LFi,-i I LF (,,2) for all program stars. 
The filled diamonds with error bars indicate detected features, 
arrows indicate upper limits. The full line represents the line 
fluxes of MWC 1080 under variable extinction. 



molecules. Differences in PAH grain sizes will change the 
ratio of number of CC bonds over number of CH bonds. 
Larger grains will have a large CC/CH ratio, which may 
lead to e.g. a smaller LF3.3/LF6.2 ratio. 

2. Dififerences in PAH ionisation can have a strong effect on 
the intensities of the PAH emission features. The PAH 
charge is ruled by the photo-ionization rate and the elec- 
tron recombination rate, and hence depends strongly on lo- 
cal physical conditions, in particular the electron density 
and the strength of the UV radiation field. 

3. Chemical differences can also affect the PAH emission 
spectrum. PAH molecules can contain other elements than 
C and H. This pollution influences the strength and espe- 
cially the peak position of the feature. The expected shifts 
in peak position due to this effect are nonetheless too small 
to be measurable in the available spectra. No definite con- 
clusions on the presence of the latter can be drawn. 

4. Dehydrogenation of PAHs, when the molecules are stripped 
of their H atoms, is believed not to play a crucial role 
JHonv et alJl200ll) . Nevertheless this phenomenon could 
add to the observed differences in the spectra. 

The 6.2 micron feature is the most frequently detected PAH 
feature in our analysis (25 times out of 46 spectra), followed by 
the 7.7 micron feature (19/46). This was to be expected, since 
these two features are both linked to the CC bond, which is an 
indispensable part of the PAH chemistry. The PAH 3.3 and 8.6 
micron features, both linked to the CH bond, are less frequent 
(12/45 and 16/46 respectively). The mean peak wavelength, 
FWHM and peak-over-continuum flux ratio of the examined 
features is given in Table The mean LF ratios of all features 
are summarized in Table E] 
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Table 7. The mean peak wavelength Ao, FWHM and peak- 
over-continuum flux ratio of the detected features. 



Feature 


<^> 


<FWHM> 


<PF/CF> 






[jum] 


Lum] 




O 4 


PAH 3.3 micron 


3.293 


0.045 


1.172 


^ 


NAN 3.4 micron 


3.415 


0.067 


1.284 


NAN 3.5 micron 


3.530 


0.059 


1.542 


b 2 


PAH 6.2 micron 


6.261 


0.198 


1.690 




PAH 7.7 micron 


7.837 


0.551 


1.575 


1 


PAH 8.6 micron 


8.641 


0.311 


1.310 


1^ 


COMP 1 1 micron 


11.254 


0.326 


1.329 


en 


Si 9.7 micron 


9.959 


2.775 


1.909 


in 



,^+ ; 


" " .. 'i +* : 




- . ^.n 


: . *"•*■•'*" • : 


'_ >- '_ 



3.2. The nanodiamond features 



0.010 0.100 1.000 

PAH 6.2 LF/CF [micron] 



o 




Only a small minority of stars in this sample display nanodi- ^ig- 1^- The LF/CF ratios of the amorphous 10 micron silicate 
amond features in their spectrum. Both examined features (at feature (Si 9.7) versus the PAH 6.2 micron feature. Diamonds 
3.4 and 3.5 /zm) were detected in the specti-a of EH as 3-1 and ^i* error bars indicate detected features, arrows represent up- 
HD 97048 (as thoroughly described bylVan Kerck hoven et alJ P^r limits. Sources with the 10 micron feature in absorption 
I2002 '). but also in the spectrum of BD-h40°4124. The spectrum h^^^ ^ negative LF/CF ratio, 
of HD 100546 possibly contains the 3.4 micron feature, while 
the 3.5 micr on feature was dete cted in the spectrum of MWC 
297 (see also'Teradaet al."200l'). 

TablesQand|^contain general information about the NAN 
features. Nevertheless, because of the small number of sources 
in this sample that display 3.4 or 3.5 micron emission, no strin- 
gent conclusions on the conditions needed to see these features 
can be drawn. We are currently carrying out a ground-based 
3.2-3.6 yum spectroscopic survey of a much larger sample of 
HAEBEs, which should allow us to draw more definitive con- 
clusions on possible correlations of the NAN and the PAH fea- 
tures in the near future. 

3.3. The amorphous 10 micron silicate feature 

In the sample of HAEBEs presented in this paper, most of the 
sources (52%) have the amorphous 10 micron silicate feature 
in emission; 6 out of 46 objects (13%) have the feature in ab- 
sorption, while it is undetected in the other cases. We note that 
at least in some cases, the presence of a 10 fim feature may be 
masked by the presence of strong PAH emission in the 8.6 and 
11.2 micron bands, so the true fraction of sources with amor- 
phous silicates is likely higher than reported here. 

The LF/CF correlation coefficients of the 10 micron fea- 
ture and the PAH features are low. In Fig. ^| the LF/CF ra- 
tio of the 10 micron feature is plotted versus that of the PAH 
6.2 micron feature. The two ratios correlate poorly, especially 
when one takes into account the upper limits (arrows) on the 
left-hand side in the plot. These arrows represent spectra in 
which the 10 micron feature is present, but where the PAH 6.2 
micron feature is undetected. Similar plots in which the other 
PAH features are plotted against the 10 micron feature show 
the same result: the strength of the emission of PAHs is uncor- 
rected with the strength of the amorphous 10 micron silicate 
feature. 

We now consider the peak-over-continuum flux ratio 
(PF/CF) of the 10 micron feature in our sample. We plotted 



PF/CF 

Fig. 11. The cumulative histogram of the PF/CF for the amor- 
phous 10 micron silicate feature. The 24 sample stars in which 
this feature is in emission are included in this plot. The lower 
filled part of the bars represents group I sources, the upper part 
group II sources. 



this dimensionless quantity for all detected 10 micron emission 
features in the sample in Fig.^] The histogram peaks around 
the average value 1.91, which means that the peak of the 10 
micron emission is on average twice as high as the underlying 
continuum. The minimum and maximum value for the ratio in 
the sample is 1.21 and 2.69 respectively. There does not seem 
to be a difference between the 10 micron emission feature in 
group I and group II sources. The average PF/CF ratios for 
the two groups are 1.85 and 2.03 respectively. Strong as well 
as more modest amorphous 10 micron silicate features appear 
in both group I and group II sources. The feature can also be 
absent in both groups. 
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Table 8. The mean line-flux ratios of the examined features. 
Indicated in the table is the mean LFa /LFb of all sample stars 
for which both A and B were detected. 



B 

A 


z 
< 

z 


in 

z 
< 
z 


X 




od 

X 


Oh 

s 

8 




PAH 3.3 


0.640 


0.310 


0.171 


0.074 


0.271 


0.138 


0.010 


NAN 3.4 




0.662 


0.119 


0.060 


0.246 


0.044 


0.004 


NAN 3.5 






0.241 


0.143 


0.416 


0.269 


0.051 


PAH 6.2 








0.437 


1.430 


0.802 


0.051 


PAH 7.7 










3.302 


1.830 


0.113 


PAH 8.6 












0.443 


0.026 


COMPll 














0.043 



3.4. The 1 1 micron feature 

The 11 micron feature is detected in 48% of the spectra. It is 
unclear how many of these detected features include the PAH 
11.2 micron band, since this feature can be blended with the 
crystalline silicate feature at 11.3 micron. Fig.^jis similar to 
Fig-El The logarithms of the ratios LFi i /LF(,,2 and LF^^/LFd^ 
are plotted. Again, the straight line indicates the variable ex- 
tinction line, based on the line fluxes of MWC 1080. Four 
sources (AB Aur, Elias 3-1, HD 100546 and HD 179218) lie 
significantly to the right of this line. This suggests that these 
sources have an 11 micron feature that contains both PAH 
and crystalline silicate emission. The presence of crystalline 
siUcates in HD 100546 and HD 179218 was already ^well- 
established from longer- wavelength ISO data ( Malfa it et alJ 
Il998l I1999I) . The detection of crystaUine silicates in AB Aur 
and Elias 3-1 is a new result from our quantitative analysis. We 
stress that the detection of a stronger than expected 1 1 micron 
feature in the four sources mentioned above does not mean that 
other sources cannot have crystalline silicates; it could just be 
masked by the presence of strong PAH emission at 11.2 fj.ni, or 
the crystalline silicates could be cooler than a few hundred K, 
causing only very weak emission in our diagnostic band at 1 1 
microns. 

When plotting the FWHM of the feature against the ratio 
log{LFi 1 /LF(,2) (Fig.ll3>. the same four sources appear to have 
much higher FWHMs. The horizontal line in the plot indicates 
the average FWHM (with exclusion of the 4 objects mentioned 
before). This also points to the presence of 2 blended emission 
features. Even so, the peak position of the 1 1 micron feature 
in AB Aur, Elias 3-1, HD 100546 and HD 179218 is compara- 
ble with the average peak wavelength (see Fig. ll4t . This might 
indicate that the peak position of the crystalline silicate feature 
varies, which can be attributed to a variety in its mineorological 
composition. We stress however that the uncertainties in these 
determinations are large. 

Sample stars that do not display significant PAH emission 
at 3.3, 6.2, 7.7 or 8.6 micron are not expected to have significant 
PAH 11.2 micron emission. Hence a detected 11 micron fea- 
ture in the spectra of these sources is likely to be solely due to 
crystalline silicate emission. In the present sample HD 142527 
(group I), HD 142666, HD 144432, HR 5999, HD 150193, HD 
163296 and VV Ser (all group II) meet these requirements. If 
these objects indeed display the crystalline silicate band in their 




log(LF_l 1/LF_6.2) 

Fig. 12. Same type of plot as Fig.|5] The logarithm of the ratio 
of the line fluxes of the PAH features at 3.3 and 6.2 micron 
\Qg{LF-},T,ILF^2) versus the logarithm of the ratio of the line 
fluxes of the features at 11 and 6.2 micron log{LFii/LF(,.2)- 
Again, the full line represents the line fluxes of MWC 1080 
under variable interstellar extinction. 
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log(LF_l 1/LF_6.2) 

Fig. 13. The FWHM of the 1 1 micron feature versus the loga- 
rithm of the ratio of the line fluxes of the features at 1 1 and 6.2 
micron log(LFi,j,/LF^2)- The dashed line represents the aver- 
age FWHM (0.214 micron) when excluding the 4 outliers. 



spectrum, then the number of sources with crystalline silicates 
in this sample is at least 1 1 (24%). 

The 1 1 micron complex does not correlate with the PAH 
6.2 micron feature, and reaches better agreement with other 
CH bond related PAHs (3.3, 7.7 micron), which is expected, 
since the 11.2 micron feature is due to a CH bond itself. The 
correlation with the PAH 8.6 micron feature, on the other hand, 
is weak. 

In Fig. ^] the ratio LF/CF of the amorphous silicate fea- 
ture versus the 1 1 micron complex is plotted. There does not 
seem to be a correlation between the two features. This can 
be interpreted in two ways. The 1 1 micron feature could be 
mostly dominated by PAH emission in this sample of stars, 
which leads to a poor correlation with the 10 micron feature 
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Fig. 14. The peak position of the detected 1 1 micron features 
versus the ratio of the line fluxes of the 1 1 micron feature over 
the 6.2 PAH feature log{LFi i /LFe.i)- The dashed line at 11.233 
micron indicates the mean peak position of the 1 1 micron fea- 
ture. 
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Fig. 15. The LF/CF ratios of the amorphous 10 micron sili- 
cate feature (Si 9.7) versus the 1 1 micron feature. Diamonds 
with error bars indicate detected features, arrows represent up- 
per limits. Sources with the 10 micron feature in absorption 
have a negative LF/CF ratio. 

like for the other PAH features (see Sect. 13. 3> . The other pos- 
sibility is that the crystalline 1 1.3 micron silicate band — when 
present in the 1 1 micron complex — does not correlate with the 
amorphous 10 micron silicate feature. Note that also the four 
'crystalline silicate stars' mentioned above do not display an 
obvious trent. This could be interpreted in terms of the second 
explanation. 

Except for Z CMa, no crystalline silicate absorption or 
emission is observed in sources which display amorphous sili- 
cate absorption (group III). The spectrum of Z CMa does dis- 
play a local minimum around 10.7-11.5 micron, which might 
be due to crystalline silicate absorption. In Fig. ^] the spec- 
trum of Z CMa is plotted. For reasons of clearness it was in- 
verted, making the absorption features appear as emission fea- 




lambda 



Fig. 16. The inverted spectrum of Z CMa, and the emission 
band spectra of HD 144432 and HD 163296. For Z CMa, 
the figure shows the inverted continuum-divided, continuum- 
subtracted flux [fconti'^) - 'F(A)\/'J^cont(A)- The other two spec- 
tra display the continuum-subtracted flux. All three spectra are 
normalized to a peak flux of 1 . An additional off'set was applied 
for the sake of clarity. The dashed line represents A- 11.3 jum 

tures. Two sources (HD 144432 and HD 163296) fliat are likely 
to display a crystalline silicate emission feature at 11.3 micron, 
are also plotted as a reference. 

4. Correlations between solid-state features and 
disk properties 

4.1. TheSED 

Using the luminosities and effective temperatures computed 
in Sect. 12.61 we plot the stars in our sample in a traditional 
Hertzsprung-Russel (HR) diagram. There is no strong correla- 
tion between the classification of the sources and the position 
in the HR diagram. In Fig. [n] the plotting symbols are pro- 
portional to the observed PAH luminosity (Lpah; Appendix). 
In this plot, theoretical pre-main-sequence (PMS) evolutionary 
tracks of st ars with masses M = 1.5, 2, 3 and 5 Mq are plotted 
(iBernasconi 1996) . Furthermore, the zero-age main sequence 
(ZAMS) is indicated by a dashed line. Most of the stars lie 
close to the main sequence (MS), which is represented by the 
dotted line. Based on this diagram, the masses of the bulk of the 
sample stars lie between 1.5 and 3 Mq. The strength of the PAH 
bands does not seem to be strongly correlated with stellar mass 
or age. Fig. ^] is a plot similar to Fig. E] for the amorphous 
10 micron silicate emission feature. The plotting symbols are 
proportional to the peak-over-continuum flux PF/LF of the fea- 
ture. The strength of this solid-state band appears to be also 
uncorrected with stellar mass and age. 

One of the main characteristics of HAEBE stars is the 
IR excess due to thermal emission of circumstellar dust. In 
Table |6l the ratio of this excess luminosity Lgxc over the stel- 
lar luminosity L is included. Ten of the sample sources have an 
IR-excess luminosity that is larger than the stellar luminosity. 
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Fig. 17. The Hertzsprung-Russel diagram for the sample stars. 
The plotting symbols are proportional to the PAH luminosity. 
Filled diamonds indicate detected emission, open symbols rep- 
resent upper limits. The full lines indicate the PMS evolution- 
ary tracks of stars with masses M = 1.5, 2, 3 and 5 Mq from 
bottom to top; the dashed and dotted lines represent respec- 
tively the ZAMS and the MS (Bernasconi 1996) . 




log(Teff [K]) 



Fig. 18. Similar plot as Fig.^] The plotting symbols are pro- 
portional to the peak-over-continuum flux of the amorphous 10 
micron silicate emission feature. 



This means that these sources cannot have passive reprocessing 
circumstellar disks. There does not seem to be a link between 
the classification of the sample stars and their Lexc 11^ ratio. The 
average value for the sources for which the ratio is smaller than 
unity, is 0.37 + 0.21 for group I, 0.39 ± 0.20 for group II and 
0.09 ± 0.04 for group III. Half of the group III objects have a 
ratio bigger than one, though, which is consistent with the idea 
that these sources are still actively accreting systems. 

In Fig. [^ the logarithm of the IR-excess luminosity 
log(Lexc) is plotted versus the logarithm of the absorbed lumi- 
nosity log(Labs)- The dashed line represents Lexc - Labs- Most 
of the sample stars lie close to this line. Exceptions include 
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Fig. 19. The logarithm of the IR-excess luminosity log(Lexc) 
versus the logarithm of the absorbed luminosity log(Labs)- 
Filled diamonds represent group I, open diamonds group II and 
squares group III sources. The dashed line is the first bisector. 
In the left upper corner, the typical error bars are represented 
by a cross. 



the group II source in the upper left corner which represents R 
CrA and the open diamond at the bottom of the plot, represent- 
ing HD 141569. The IR excess luminosity of R CrA is due to 
confu sion with the nearby source CrA IRS7 (Choi & TatematsU; 
I2OO4I) . Furthermore, a set of group I sources with low Labs val- 
ues, display relatively high Lexc values. An explanation for this 
extra IR excess luminosity could be the presence of a late-type 
companion. Also the orientation of the system can play a role, 
especially in group I sources, where the flaring of the disk im- 
plies that the line of sight to the central star already passes 
through the disk at fairly low inclinations (i ~ 45°). 

4.2. The SED versus the spectral features 

The different parameters that characterize the SED were com- 
pared to the parameters that describe the infrared solid state 
bands. In this section we summarize the results. 

The UV luminosity of the central star correlates with 
the PAH luminosities: the absolute PAH luminosity increases 
with increasing UV radiation. The ratio Lpah/^uv however de- 
creases with increasing UV luminosity. In Fig. |20l the latter 
correlation is demonstrated. The dashed line in this figure rep- 
resents log(LpAH/Luv) - -0.395 log(Luv) - 1-99 or equivalently 
LpAH = 0.01 L^f. If the emitted PAH luminosity would in- 
crease linearly with increasing UV radiation of the central star 
(log(LpAH '^ i'uv), the ratio Lpah/Li„ would have been con- 
stant for all Lyy. This is obviously not the case. A higher stellar 
UV luminosity indeed increases the PAH luminosity, but does 
this in a sub-linear way: the PAH luminosity increases with 
Luv where p = 0.6 < 1. It appears that either the efficiency 
of the absorption/emission process decreases with increasing 
UV strength, or that the increasing hardness of the UV photons 
plays a role. Since hotter stars emit not only more UV photons, 
but also relatively more high-energy UV photons, the observed 
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Fig. 20. The PAH-over-UV luminosity log(LpAH/iuv) versus 
the intensity of the UV radiation field log(Luv)- The PAH lumi- 
nosity is the sum of the luminosities of the PAH features at 3.3, 
6.2, 7.7 and 8.6 micron. In the left upper corner, the cross indi- 
cates the typical errors. FiUed diamonds refer to group I, open 
diamonds to group II and squares to group III objects. Arrows 
represent upper limits. The dashed line represents the best lin- 
ear fit to the data: log(LpAH/iuv) = -0.395 log(Luv) - 1-99. 

trend might be a contrast effect. Short-wavelength UV photons 
cannot be absorbed by PAHs (no increase of Lpah), but do con- 
tribute to the total UV luminosity Luv, hence decreasing the 
luminosity ratio. 

It is quite remarkable that some group II sources, even 
though they have comparable UV luminosities as their group I 
counterparts, do not display PAH emission. This indicates that 
the appearance of PAH emission is linked to the group I/II clas- 
sification of the sources, and does not depend on L^v only. 

Suppose that the IR emission is due to PAH molecules that 
are homogeneously distributed in a spherical, optically thin 
halo around the star, which has the radius of the ISO-SWS 
beam (11") at the appropriate distance. It is then possible to 
estimate the mean particle density of the PAH molecules (p; 
see Appendix) in that halo. In Table|9lthe computed values are 
listed. The detected emission of the PAH features at 3.3, 6.2, 
7.7 and 8.6 micron has been included. Note that not all emitted 
PAH flux (e.g. the PAH 11.2 micron flux) is accounted for. This 
implies that — under the given assumptions — the computed 
molecule density is a lower limit. For most objects, the result- 
ing values are high compared to the typical ISM PAH density 
(ptyp = 5 X IQ-^ cm--'; Appendix). We interpret this as an in- 
dication that the PAH emission does not emanate from a halo 
around the central star For most stars, a denser environment is 
ne eded to explai n the observed emission. 

iHabart et alJ (12004) have modelled PAH emission in 
flared circumstellar disks. In their models, they assume PAH 
molecules containing 40-100 carbon atoms. For PAHs of these 
sizes, the influence of ionization and photo-evaporation affects 
the line-flux ratios of the PAH features, in particular the rela- 
tive strength of the 3.3 micron feature. Even though a clear de- 
pendence of the PAH line fluxes on the stellar UV luminosity 



Table 9. Supposing a homogeneous distribution of the PAH 
molecules in a halo around the central star, with a radius of the 
size of the ISO-SWS beam (11") at the appropriate distance, the 
PAH particle density p is computed based on the Lpah/^uv ra- 
tio. All objects in which PAH emission is detected are adopted 
in this table and are listed according to decreasing PAH-over- 
UV luminosity. 



Does PAH emission emanate from 


a lialo? 


Object 


Group 


-^PAh/^^uv 


P 
[cm-3] 


Elias3-1 




0.36 


5.51 10-* 


RRTau 




2.5 10-2 


3.82 10-5 


HD 34700 




1.1 10-- 


8.13 10-5 


HD 97048 




1.1 10-- 


1.51 10-5 


HD 34282 




1.0 10-2 


6.21 10-* 


VXCas 




7.7 10--3 


3.00 10-*^ 


HD 169142 




6.8 10-3 


1.14 10-5 


HD 100546 




6.4 10-3 


1.51 10-5 


HD 179218 




5.5 10-3 


5.60 10-« 


HD 100453 




5.2 10-3 


1.14 10-5 


HD 95881 




4.2 10-3 


8.70 10-^ 


Wra 15-1484 




3.5 10-3 


1.14 10-* 


HD 139614 




3.3 10-3 


5.78 10-*^ 


VVSer 




2.6 10-3 


1.93 10-'^ 


AB Aur 




2.2 10-3 


3.87 10-'^ 


HD 142666 




2.0 10-3 


3.42 10-« 


HD 135344 




1.6 10-3 


2.86 10-^ 


CD-42° 11721 




1.3 10-3 


8.16 10-^ 


HD 144432 




9.8 10-* 


1.65 10-'^ 


BD+40°4124 




9.3 10-* 


2.31 10-^ 


HD 31648 




8.5 10-* 


1.58 10-*^ 


HD 163296 




7.1 10-* 


1.42 10-*^ 


HD 142527 




6.9 10-* 


1.15 10-*^ 


HD 141569 




5.3 10-* 


1.31 10-*^ 


MWC 1080 




1.4 10-* 


1.54 10-'* 


MWC 297 


III 


6.8 10-5 


6.64 10-** 


HD 200775 




3.1 10-5 


1.72 10-** 



was observed in our present sample, no correlations between 
the relative strength of the PAH features and the UV radiation 
were noted. The line-flux ratio of the 3.3 and 6.2 micron feature 
-which ranges from 9 to 94%- is independent of the central 
star's UV field. 

The excesses in the K band, at 60, 850 and 1300 micron 
were compared to the PAH data. Fig. 1211 shows the luminos- 
ity of the PAH 7.7 feature versus the excess at 60 micron. 
The sources with the faintest 60 micron excesses (by definition 
mostly group II objects) are those that have the faintest PAH 
emission. Other PAH features display comparable behaviour. 

Notwithstanding these correlations, in general the lumi- 
nosities of the IR emission bands do not correlate well with the 
other measured excesses. The lack of correlation between the 
850 or 1300 micron excess and the strength of the PAH features 
shows that the disk mass has no influence on the resulting PAH 
spectrum. This is consistent with the idea that PAH emission 
emanates from the surface layers (atmosphere) of the circum- 
stellar disk, which contains only a small fraction of the disk 
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Fig. 21. The PAH 7.7 micron luminosity logiL^AWin) versus 
the excess at 60 micron (in magnitudes). Filled diamonds rep- 
resent group I, open diamonds group II and squares group III 
sources. Small symbols indicate undetected features, arrows 
stand for upper limits. In the left upper corner the typical er- 
ror bars are displayed. 

mass. The FWHMs of the IR PAH bands do not correlate with 
the excesses either. 

The Rayleigh-Jeans part of the SED can be modelled with 
a power law, \og(A'F;i) ~ A^" in which n is the (sub-)mm spec- 
tral index. There is no clear correlation between this index and 
the strength of the PAH features. The sub-mm slope of the 
spectrum can be a proxy for the grain size distribution of the 
cold, large grains in the outer parts and mid-plane of the disk 
JAcke et al. 2004, henceforth A04). This non-correlation indi- 
cates that the grain size distribution of the bulk of the disk's 
material has no influence on the PAH emission. This is again in 
agreement with the hypothesis that PAH molecules are excited 
in and radiate from the disk's atmosphere. 

The parameters that describe the amorphous silicate 10 mi- 
cron feature do not correlate with the IR flux excesses. This 
indicates that the strength of the 10 micron feature does not 
depend on the disk mass. The LF and FWHM of the silicate 
band do not correlate with the sub-mm spectral index either. 
The grain size distribution of the cold particles in the outer parts 
of the disk and the warm small silicate grains are independent 
of each other. Nevertheless, if vertical mixing in the disk is ef- 
ficient, one may expect that only the smallest grains make it 
to the surface. Hence the amorphous 10 micron silicate feature 
may be indicative for the smallest grain sizes (still) present in 
flie disk (A04). 

We compared the SEDs of the four stars for which the 1 1 
micron feature is a blend of the PAH 11.2 micron band and 
the crystalline silicate feature at 11.3 micron (Elias 3-1, AB 
Aur, HD 100546 and HD 179218) with the SEDs of die other 
sample stars. The sub-mm spectral index n, 10 micron PF/CF 
ratio and 1300 micron flux excess of these group I sources 
make them normal members of this group. The objects AB 
Aur, HD 100546 and HD 179218 have steep sub-mm slopes 
(n~4.20), high 10 micron PF/CF values (larger than the aver- 



age value 1.91), and intermediate 1300 micron excesses (~10 
mag). The circumstellar material of these sources thus con- 
tains cold grains in the outer parts of the disk, significantly 
smaller than mm-siz ed (AV04) an d small (~ 0.1 um) warm sili- 
cate grains (Bouwm an et al.l200lH van Boek el et al.l2003l) . The 
estimated disk mass is of the order of 1 % of the stellar mass. 
Elias 3-1 on the other hand is somewhat more peculiar, since 
n is equal to 3.22 and the 10 micron peak-over-continuum flux 
ratio is 1 .70. The circumstellar matter in this system seems to 
be more evolved. The disk mass of this source is also high com- 
pared to the stellar mass (30%). This is cl ose to the the oretical 
upper limit for gravitational stability ('e.g. lGammiel2001l) . 



5. Interpretation 

We use plotting diagrams like Fig.|5l and include the strength 
of the IR spectral features by scaling the plotting symbol pro- 
portional to LF/CF. In Fig. |22] the sum of the LF/CF ratios 
of the PAH features at 3.3, 6.2, 7.7 and 8.6 micron is indicated 
in this manner. The filled diamonds in the plot represent de- 
tected features, while open diamonds stand for upper limits. 
Squares refer to group III sources. The sizes of the diamonds 
and squares represent the strength of the PAH features on a 
logarithmic scale. We specifically indicate group III sources 
since, because of their location in the diagram, they could be 
confused with group I sources. That group III sources appear 
in the group I part of the plot is to be expected, since these 
sources are believed to be in an earlier evolutionary stage; the 
objects are stifl highly embedded in the circumstellar environ- 
ment, which blocks out the stellar light. This makes their SEDs 
very red and hence they show up in the lower right part of the 
diagram. 

It is striking that the sources that display the strongest PAH 
emission are mostly group I sources. The .Dullemond (.20021) 
model suggests that these objects have circumstellar disks with 
a flaring geometry. PAH molecules in the atmosphere of the 
flaring part of these disks are directly irradiated by the UV pho- 
tons of the central star. This implies that the molecules will be 
excited and will radiate in the IR, no matter how far out in the 
disk, as long as the UV field is strong enough. PAH molecules 
in the atmosphere of the self-shadowed, geometrically flatter 
disks of group II sources barely see any direct UV radiation, be- 
cause the puff'ed-up inner rim blocks out most of the starlight. 
Hence, in these systems the IR emission bands are expected to 
be much fainter Most objects in our sample are consistent with 
this hypothesis. 

Except for MWC 297, no PAH emission is detected in 
group III sources. The circumstellar matter of these highly 
embedded objects can probably accurately shield the stellar 
UV radiation close to the central star The volume of ex- 
cited PAH molecules will likely be too small to be observable. 
Furthermore, the PAH IR radiation has to travel through the 
dense circumstellar environment. This prohibits the appearance 
of PAH emission in the group III spectra. 

In Fig. |23lthe LF/CF ratio of the amorphous 10 micron 
silicate feature is plotted in the Lnir/^ir versus [12] - [60] di- 
agram. Again, the filled diamonds refer to the detected emis- 
sion features, while the open diamonds represent upper limits. 
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Fig. 22. Similar plot as Fig.|5] The plotting symbols are scaled 
proportional to the strength (2 LF/CF) of the PAH features 
at 3.3, 6.2, 7.7 and 8.6 micron. The outer limit of the par- 
tially filled plotting symbol represents the upper limit for the 
strength; the filled inner part refers to the part that was actu- 
ally detected. Diamonds represent group I and group II sources, 
squares group III objects. 

Group III sources, which display amorphous silicate absorp- 
tion, are not included in this plot. Note that the silicate emis- 
sion feature appears in both groups. When present, the strength 
of the emission feature seems to be independent of the classi- 
fication of the source. Nevertheless, a larger fraction of group 

I sources have an undetected amorphous silicate feature (10/18 
versus 6/22 in group II). MOl suggested that this is probably a 
selection effect; group II objects that do have 10 micron silicate 
emission are brighter, and hence are observationally favoured 
over group II sources without siUcate emission. This selection 
effect plays a smaller role for group I sources, since they are by 
definition brighter IR emitters. ^^^^_^^^^^ ^_^ 

According to the definition of IPuUemond et alJ (l2003l) . 13 
of our sample sources are UX Orionis (UXOR) stars. In their 
paper, the authors argue that the disk geometry of this class of 
stars must be self-shadowed and hence that UXORs are group 

II objects. In the present paper we have shown that PAH emis- 
sion in group II sources is weak. Of the 13 UXORs in our sam- 
ple, only 3 sources do show significant PAH emission: VX Cas, 
HD 34282 and RR Tau. They ai-e labelled in Fig.|24| We note 
that all three lie close to the empirical separation line between 
group I and II; they might be transitional objects in which the 
shadow cast by the inner rim is relatively small. We remind the 
reader that R CrA and LkHa 224 are not in the plot, and that 
these sources were classified based on their UXOR behaviour. 

In Fig.|25j a selection of 14 of our sample stars shows that 
the 10 micron emission feature tends to become broader when 
the peak-over-continuum flux ratio decreases. SmaU warm sili- 
cate grains (~0. 1 jim) cause a distinct feature that peaks around 
9.7 micron, while larger grains (~2 jum) induce a broader, less 
pronounced feature wi t h a peak towards longe r wavelengths 
(Bouwman et al. 20011) . Ivan Boekel et al.l 12003) interpret this 
sequence in terms of grain growth in the disk. In Fig. 1251 this 



: Group 1 

> 


* ♦ ♦ 

; 
I 


♦ 

♦ 


o - 
♦ 










Group 1 



[12]-[60] 

Fig. 23. Similar plot as Fig.|22lfor the amorphous 10 micron 
silicate feature. The plotting symbols are scaled proportional 
to the LF/CF ratio of the feature. The filled diamonds indi- 
cate stars with a detected 10 micron emission feature, the open 
diamonds indicate upper limits. 
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Fig. 24. Similar plot as Fig [S] The filled symbols indicate 
UXORs, the open symbols are the other HAEBEs in this sam- 
ple. The 3 UXOR stars that display clear PAH emission are VX 
Cas, HD 34282 and RR Tau. 



sequence is reproduced. Our analysis confirms observationally 
the correlation between the shape and the strengt h of the amor- 
phous 10 micron silicate band. In their Letter, van Boekel et alJ 
1I2OO3I) only plotted group II sources. In this study we have 
shown that the 10 micron feature is independent of the clas- 
sification of the objects. Therefore, we also allowed group I 
sources in Fig.|25] The correlation plot (Fig.l26ll of the FWHM 
versus the PF/CF ratio of the amorphous 10 micron silicate 
feature shows that indeed the FWHM is larger when the peak- 
over-continuum flux ratio is lower Again, no difference is ob- 
served between group I and group II sources. The dashed line 
represents the best fit and has a slope of -0.38 micron. 
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Fig. 25. The amorphous 10 micron siHcate band for 14 of the 
sample stars. The figure shows the continuum-subtracted flux 
'F(A) - ^FcontC'^), which was divided by the average value of the 
underlying continuum {!Fcont('t))- The spectra of the stars are 
sorted by increasing peak-over-continuum flux from a to n: a-R 
CrA; b-HR 5999; C-V376 Cas; d-HD 200775; e-T CrA; f-VX 
Cas; g-HD 104237; h-BF Ori; i-HD 31648; j-HD 163296; k- 
SV Cep; 1-HD 144432; m-UX Ori; n-AB Aur. For clarity, the 
zero levels of the 14 spectra are shifted. 
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Van Boekel et al. f2003j suggested the amorphous 10 mi- 
cron silicate emission originates from the disk's surface. It 



PF/CF 

Fig. 26. The FWHM of the amorphous 10 micron silicate fea- 
ture versus the peak-over-continuum flux ratio PF/CF. Filled 
and open diamonds refer to group 1 and group 11 sources re- 
spectively. The dashed line represents the best fit to the data; 
FWHM = -0.38 X PF/CF -i- 3.67. The typical error bar is indi- 
cated in the upper left corner 



might appear contradictory that the strength of PAH emission 
does depend on the shape of the SED (which represents the disk 
geometry), while the 10 micron feature does not. Nevertheless, 
the excitation mechanisms for the two emission sources are 
not the same: the non-equilibrium PAH emission occurs during 
temperature fluctuations after absorption of a UV photon while 
the larger silicate grains are in thermal equilibri um with the 
radiation field of the central star. The models of Habart et alJ 
( 12004) show that the PAH emission emanates mostly from the 
outer parts of the disk (~100 AU), while the warm silicate 
emission is confined to the innermost disk regions (~ a few 
AU). This reflects the diff'erent excitation mechanisms of the 
small PAH molecules and larger silicate grains, and is in accor- 
dance with the observations: spatially resolved PAH emission, 
on scales of ~10-100 AU , has been detected arou nd HAEBEs 
bv lGeersetaD(l2004 and lvan Boekel et al.l ( l2004 . 

The results of the present paper on both the PAH emis- 
sion and the silicate emission are consistent with a scenario 
in which Group I sources have disks that are flaring, whereas 
Group II sources have flatter disks that are shadowed by the 
disk's pufFed-up inner rim. 



6. Conclusions and discussion 

We summarize the conclusions which were drawn in the previ- 
ous sections. 

1. Although in general the PAH emission-line spectrum of 
HAEBEs is fairly homogeneous, we find clear evidence of 
differences in size, chemistry and/or ionization for the PAH 
molecules within the group of Herbig Ae/Be stars. 

2. The 6.2 micron feature is the most frequently detected PAH 
feature in our analysis. However, in 56% of the spectra 
where the 6.2 micron feature is detected, the PAH 3.3 mi- 
cron feature is Mndetected. Hence the 3.3 micron feature 
is not a good indicator for the presence of PAH emission, 
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since quite a large fraction of the PAH emitters may be 
overlooked if one only considers this feature. 

3. PAH emission occurs more frequently, and is stronger in 
MOl Group I sources than in Group II sources. If we as- 
sume that PAH emission originates in the disk's atmo- 
sphere, this result provides strong support for the interpre- 
tation of Group I sources as sources with flared disks, and 
Group II sources as those with self-shadowed disks. 

4. The strength of the PAH bands is correlated with the 
strength of the UV radiation of the central star, but the 
relative strength of the PAH emission decreases with in- 
creasing Luv This might be due to the fact that the power 
absorbed by PAHs changes with the spectral energy distri- 
bution of the stellar radiation field and decreases with in- 
creasing hardness. 

5. Except for MWC 297, no PAH emission is observed in 
group III sources. 

6. Nanodiamond features were detected in four sample stars, 
of which three are group I objects and one is group III. 

7. The shape and strength of the 10 micron feature is inde- 
pendent of the MOl classification; strong, weak and even 
absent emission features are observed in both groups. The 
10 micron feature in the spectra of the sample sources 
indicates that warm amorphous grains are processed in 
HAEBEs. 

8. Neither the PAH features nor the 10 micron feature depend 
on disk mass. The masses of the disks are also independent 
of the classification. 

9. The 11 micron feature does not correlate with the amor- 
phous silicate feature, which could indicate that the PAHs 
are the dominant ingredient of the feature and/or that crys- 
talline siUcates have an independent evolution. The 1 1 mi- 
cron complex is in general stronger in group I sources. 

10. Crystalline silicates in emission occur in both types of 
sources. Z CMa is the only source where the feature is pos- 
sibly seen in absorption. 

We confirm in the present paper that the strongest PAH- 
emission occurs in group I sources, as was first noted by MOl. 
This is consistent with the idea that group I sources have a 
flared disk geometry, while group II so urces are systems with 
self-shadowed disks ( Dulle mondl2002l) . In a flared disk geom- 
etry, the UV photons of the central star can reach the PAH 
molecules in the surface layers of the disk in a direct way. PAH 
molecules in the atmosphere of self-shadowed disks are shaded 
from the UV field by the puff'ed-up inner rim which casts its 
shadow over the outer parts of the disk. Therefore, the PAH 
spectrum of group II sources is expected to be much fainter 
than that of group I sources. Most sample stars follow this hy- 
pothesis. 

The amorphous 10 micron silicate feature in this sample 
of HAEBE stars displays the characteristics of grain growth, 
independent of the classification of the sources. The small hot 
silicate grains from which this feature originates seem to have 
little interaction with the PAH molecules, even though they are 
believed to radiate from the same locus: the disk atmosphere. 
Crystallinity seems to be present in more than 1/4 of the stars 
in our sample. 



The global interpretation of the SED of HAEBEs in terms 
of disk geometry is in many ways consistent with the obser- 
vations. Future work with new generation IR instruments (e.g. 
Spitzer) and interferometers (e.g. MIDI/ AMBER on the VLTI) 
will provide the community with new high-quality observa- 
tions of HAEBE stars, which will offer a good test for the hy- 
potheses suggested in this paper 

Appendix 

We estimate the density of PAH molecules in a spherical, op- 
tically thin halo of radius R around the central star, supposing 
a homogeneous distribution of the PAHs. Furthermore, we as- 
sume that all UV radiation absorbed by the PAH molecules is 
re-emitted in the IR bands. 

The cross-section of a single PAH r nolecule with 50 C 
atoms is typicafly o-c = 2.5 x 10 ''^ cm^ (Li & Drainell2002i) . 
The luminosity of the PAH emission Lpah depends on the num- 
ber of emitting PAHs A^pah, the average energy of the absorped 
UV photons (/iuv), the number of UV photons A^uv that pass by 
per second and the cross-section of a single PAH molecule over 
the average surface S of the sphere: 

4-7T , 

= — /?Vx<-Euv)x 



O-C 



<£uv) S 



^ = 0.37p[cm-^^]7;[AU] 



(1) 



Assuming that all IR PAH emission is observed in the mea- 
sured PAH bands, the sum of the line fluxes of these features is 
connected to the PAH luminosity by 



LpAulLe] = 2.4 10*' ^ Anidlpcjf X LF j 

PAH j 
PAH j 



w_ 



(2) 
(3) 



in which the terms of the sum are the luminosities of the indi- 
vidual PAH features. Hence from formulas[nand|2l an estimate 
for p can be derived, assuming that the radius of the halo R is 
equal to the radius of the ISO-SWS beam (11") at the distance 
d of the source, or 

R[Al]] = lldlpc] 

A typical value for the ISM abundance of PAHs over H 
atoms ([PAH/H] = 10"^) can be found in Tielens et alJlll999l) . 
Typical ISM electron densities pg- are in the order of 5 cm ■'. 
Hence a typical value for the ISM density of PAH molecules is 



Ptyp - 



PAH 



H 



xpe- = 5 X 10 ^cm 



(4) 
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